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1.  INTRODUCTION 

A number  of  time-mean  velocity  profile  models  have  been  proposed 
for  the  two-dimensional,  constant  property,  turbulent  boundary  layer. 

The  purpose  of  this  report  is  to  present  comparisons  of  several  such 
models  with  experimental  data  and  to  offer  a critical  evaluation  of  each 
model  based  on  these  comparisons.  It  should  be  stated  at  the  outset 
that  while  poor  comparisons  with  data  constitute  a strong  reason  for 
rejecting  or  modifying  a given  model,  no  amount  of  good  data  fits  can 
ever  prove  the  validity  of  that  model.  However,  favorable  agreement 
with  data  over  a wide  range  of  flow  conditions  gives  a degree  of  confi- 
dence in  a given  theory  and  seems  the  logical  first  step  in  the  develop- 
ment of  a reliable  prediction  method.  Consequently  the  main  objective 
here  is  to  observe  how  well  several  velocity  profile  models  compare  with 
experimental  data. 

Specifically  four  models  will  be  considered.  The  Van  Driest  (1956) 
model  (with  and  without  its  various  modifications  for  pressure  gradients) 
and  the  unsteady  wall  layer  model  (Walker  et  al.,  1976;  Walker  & Abbott, 
1977;  Walker  & Scharnhorst,  1977)  will  be  compared  in  the  inner  region 
in  §3.  In  the  outer  region  Coles'  (1956)  "law  of  the  wake"  and  an  outer 
region  similarity  profile  (discussed  by  Mellor  and  Gibson  (1966)  and 
Fendell  (1972))  will  be  considered  in  §4  and  §5  respectively,  each  in 
conjunction  with  the  unsteady  wall  layer  inner  model.  A detailed 
description  of  each  model  will  be  delayed  until  the  appropriate  section. 


2 


All  four  models  apply  to  a two-dimensional,  constant  property, 
turbulent  boundary  layer  which  is  steady  in  the  time-mean  sense.  The 
governing  equations  are  mass  continuity  and  the  streamwise  momentum 
equation  given  by. 


3u* 

3x* 
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u*  1^  + V* 
3X* 


3u^ 

3y^ 
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= U 


e dx’' 


+ V 


3^u* 


3y 


★ 2 


3y^ 


(u'v* ), 


(1.1) 

(1.2) 


Here  (x*,  y*)  are  Cartesian  coordinates  with  corresponding  time-mean 
velocities  (u*,  v*)  and  fluctuating  components  (u',  v');  v is  the  kine- 
matic viscosity.  The  Reynolds  normal  stresses  in  (1.2)  are  neglected 
thereby  excluding  relaminarization  or  incipient  separation  where  the 
gradient  with  respect  tox*of  these  terms  may  possibly  become  important. 

The  Reynolds  shear  stress  term,  -u'v',  marks  the  distinguishing  feature 
between  the  laminar  and  turbulent  boundary-layer  equations.  The  func- 
tional form  of  both  the  normal  and  shear  stress  terms  is  unknown  and  it 
is  this  Information  which  is  lost  in  the  time  averaging  process.  It  has 
been  traditional  in  the  past  to  attempt  to  postulate  the  form  of  these 
unknown  functions  based  on  the  belief  that  constitutive  relations  exist  re- 
lating the  Reynolds  stresses  to  the  mean  profile.  Alternative  approaches, 
such  as  turbulent  kinetic  energy,  obtain  equations  for  the  Reynolds 
stress  terms  but  at  the  considerable  expense  of  introducing  additional 
unknown  functions. 

Before  moving  to  the  data  comparisons  there  are  two  points  which 
deserve  attention.  The  first  concerns  a feature  common  to  all  four 
models  mentioned  above.  Each  contain  independent  parameters  which  can 
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be  adjusted  to  alter  the  profile  form  and  one  objective  of  the  data 
comparisons  is  to  determine  what  set  of  parameter  values  best  fits  the 
given  data  and  to  discern,  if  possible,  any  trends.  However  there  is 
one  quantity,  namely  the  time-averaged  friction  or  shear  velocity,  u^, 
which  appears  either  explicitly  or  implicitly  in  all  the  profile  models 
and  which  will  be  treated  here  as  a parameter.  There  may  be  some  ques- 
tion as  to  why  this  is  considered  a parameter  since  experimental  values 
are  always  given.  The  primary  reason  concerns  the  method  of  determining 
u^  from  experimental  data.  An  excellent  review  of  various  techniques  is 
given  by  Brown  and  Joubert  (1969)  and  a critical  discussion  of  those 
techniques  and  others  can  be  found  in  Walker,  et  al.  (1976).  Basically 
the  methods  fall  into  three  classifications:  graphical,  indirect  and 

direct  measurements. 

The  graphical  techniques  include  the  classical  Clauser  crossplot 
method  and  a secant  method  described  by  Kline  et  al.  (1967).  The 
Clauser  crossplot  method  assumes  the  existence  of  a universal  "law  of 
the  wall"  or  "log-law": 

= ^lny^  + C.  , (1.3) 

where 


This  relation  is  used  in  conjunction  with  the  measured  data  to  plot  i7* 
versus  y*.  A determination  of  u^  is  made  from  the  intersection  of  the 
experimental  profile  and  the  relation  given  by  (1.3).  To  implement  this 
technique  k and  are  usually  assumed  constant,  with  0.41  and  5.0, 
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respectively,  being  typical  values.  However,  suggested  values  for  k 
vary  from  0.39  by  Schlichting  (1968)  to  0.44  by  Kays  et  al.  (1970). 
Suggested  values  range  from  4.9  by  Cebeci  and  Smith  (1974),  and 
Mellor  and  Gibson  (1966)  to  5.56  by  Schlichting  (1968).  Empirical 
relations  for  k and  as  functions  of  Reynolds  number  have  even  been 
proposed  by  Simpson  (1970).  Clearly  such  variations  are  not  palatable 
if  accurate  and  reliable  values  for  are  to  be  obtained.  In  general, 
C^.  for  a boundary  layer  will  depend  on  the  flow  conditions  and  should 
be  expected  to  be  a function  of  x*  (Fendell,  1972);  the  assumption  that 
C.  is  a universal  constant  is  not  supported  by  Walker  et  al.  (1976)  or 
the  data  comparisons  in  this  report. 

The  secant  method  uses  the  slope  of  a line  through  data  points 
near  the  wall  to  estimate  the  slope  at  the  wall.  A value  of  u^  is  then 
calculated  from  its  definition; 


This  technique  is  questionable  since  near-wall  data  may  not  be  reliable 
due  to  probe  interference;  moreover  this  technique  is  an  extrapolation 
procedure.  Extrapolation  is  always  a dangerous  procedure  but  would  seem 
especially  so  here  in  view  of  the  rapid  variation  of  the  profile  near 
the  wal 1 . 

Preston  tubes  and  flush-mounted  surface  hot  film  probes  represent 
indirect  techniques  for  determining  u^.  However,  they  rely  on  the  law 
of  the  wall  for  calibration  and  specifically  on  the  assumption  that  < 
and  are  universal  constants.  Thus  they  have  the  same  deficiencies  as 
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Clauser's  crossplot  method  and  of  course  give  nearly  the  same  results 
(Strickland  & Simpson,  1973). 

The  only  direct  measurement  instrument  is  the  floating  element  as 
described  by  Brown  and  Joubert  (1969).  Even  though  it  appears  to  offer 
promise,  it  is  subject  to  secondary  pressure  forces  due  to  the  gap 
around  the  element  and  flow  through  its  cavity.  This  makes  calibration 
very  difficult  but  it  has  been  attempted  by  Brown  and  Joubert  (1969), 
Pierce  and  Zimmerman  (1972)  and  Miller  (1972). 

In  view  of  these  uncertainties  in  determining  u^,  it  will  be  con- 
sidered an  independent  parameter  which  can  be  adjusted  to  generate  the 
best  fit  to  a given  data  set.  At  the  same  time,  although  the  quoted 
experimental  values  are  considered  to  be  rough  estimates  of  u^,  large 
discrepancies  between  the  experimental  and  optimized  value  obtained  in 
this  study  would  clearly  be  unacceptable.  The  tabulated  data  (to  be 
given  subsequently  in  this  report)  reveal  that  the  differences  in  the 
two  u^  values  are  not  large  but  are  significant. 

It  remains  to  define  the  basis  for  determining  a best  fit  and  the 
best-fit  criterion  is  the  second  point  which  deserves  attention  here. 
For  this  study  a root-mean-square  error  was  defined  as: 


1 ^ 
l 


N ^ 

^ n=l 


‘^Data^'^n^  ~ ^Analytical 


where 


N E number  of  data  points, 
u^ata  ^ experimentally  measured  velocity  at  y*, 
Analytical  ^ appropriate  analytical  profile  at  y*, 

Ug  = freestream  velocity. 


(1.5) 
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An  error  criterion  based  on  the  normal  distance  between  the  data  point 
and  the  analytical  profile  could  be  defined  and  indeed  was  attempted  by 
York  and  Abbott  (1973).  However  there  is  an  iteration  involved  which 
was  not  implemented  by  York  and  Abbott  (1973)  and  this  resulted  in  an 
unnatural  weighting  being  assigned  to  the  near  wall  points.  It  is  felt 
that  the  standard  error  criterion  given  by  (1.5)  represents  an  objective 
basis  of  comparison  and  in  addition  it  parallels  the  technique  used  in 
the  Stanford  Conference  (Coles  & Hirst,  1969).  The  best  fit  to  a given 
experimental  profile  is  defined  as  that  set  of  profile  parameter  values 
which  minimizes  e as  defined  by  (1.5). 

The  plan  of  the  report  is  as  follows.  The  particular  optimization 
techniques  will  be  discussed  in  Section  2.  The  comparison  of  the  two 
inner  region  profiles  will  be  presented  in  Section  3.  In  Sections  ^ 
and  5 the  two  outer  region  models  will  be  tested  and  a critical  evalua- 
tion of  each  model  will  conclude  this  report  in  Section  6. 
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2.  OPTIMIZATION  TECHNIQUES 


To  minimize  e one  of  two  different  optimization  techniques  was  used 
depending  on  the  model  to  be  tested.  For  the  Van  Driest,  unsteady  wall 
layer  and  Coles'  models,  the  efficient  Newton-Greenstadt  gradient  method 
was  applied,  since  each  model  can  be  written  as  an  analytical  function. 
However,  the  similarity  profile  discussed  by  Fendell  (1972)  involves  the 
numerical  solution  of  a differential  equation  and  here  a direct  search 
technique  proved  more  feasible.  The  important  features  of  each  tech- 
nique are  summarized  below. 


2.1  Newton-Greenstadt 

The  Newton-Greenstadt  minimization  method  is  a variation  of  a 
general  iterative  technique  and  is  described  in  detail  by  Bard  (1974). 
The  Newton  nomenclature  signifies  that  a quadratic  approximation  to  the 
objective  function  (c)  is  made  at  an  estimate  of  the  minimum  at  x = x^. . 
Here  x is  a vector  which  has  as  components,  the  parameters  which  are 
allowed  to  vary  in  the  optimization  procedure.  The  necessary  condition 
for  the  existence  of  an  extremum  then  becomes: 


- «i  11  = Xj  ♦ X'i 


(2.1) 


where  for  the  i^^  iteration: 


x^.  = the  parameter  vector, 

H^  = the  Hessian  (matrix  of  second  derivatives  of  e). 
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q^.  E the  gradient  vector,  8e/9x^, 
v^.  = the  search  direction. 

Equation  (2.1)  is  often  referred  to  as  Newton's  method  and  in  practice 
to  assure  convergence  of  the  method,  it  is  necessary  for  to  be  posi- 
tive definite.  However  this  may  not  be  the  case  in  general,  except 
as  the  minimum  is  approached;  if  there  is  a small  negative  eigenvalue 
associated  with  H^.  = H(x^. ) then  the  step  given  by  (2.1)  will  be  such 
that  optimization  search  will  be  driven  away  from  the  minimum.  On  the 
other  hand  if  has  a small  positive  eigenvalue,  then  equation  (2.1) 
implies  a large  step  in  the  parameter  direction  associated  with  that 
eigenvalue;  this  situation  is  often  not  desirable.  The  modification 
of  Newton's  method  to  overcome  these  difficulties  that  was  used  here 
is  due  to  Greenstadt  and  is  discussed  by  Bard  (1974).  The  basis  of 
the  method  is  to  expand  the  Hessian  in  terms  of  its  eigenvalues  and 
eigenvectors.  If  a negative  eigenvalue  is  encountered  it  is  made 
positive  by  taking  the  absolute  value;  if  the  ratio  of  the  largest  to 
the  smallest  eigenvalue  is  too  large  (>  10®),  the  smallest  eigenvalue 
is  arbitrarily  increased.  With  these  redefinitions  of  the  eigenvalues 
of  H^.  (when  the  problem  occurs),  a revised  Hessian  may  be  calculated  in 
terms  of  the  old  eigenvectors  and  new  eigenvalues.  The  object  of  this 
modification  is  to  ensure  that  the  optimization  will  always  proceed  in 
the  right  direction  and  also  that  it  will  not  be  misled  into  taking  too 
large  a step  in  one  of  the  parameter  directions.  The  iteration  prescrip- 
tion (2,1)  can  now  be  used  safely  to  execute  the  algorithm.  At  any 
stage  in  the  iterative  procedure,  the  step  indicated  by  (2.1)  may  be 
too  large  which  can  be  rectified  by  reducing  the  magnitude  of  the 
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vector  until  a successful  step  has  been  taken.  Here  a successful 
step  is  defined  as  one  which  leads  to  a smaller  value  of  the  objective 
function.  Once  a successful  step  was  taken,  a process  of  either 
interpolation  or  extrapolation  was  used  to  accelerate  the  progress  of 
the  algorithm.  This  general  technique  is  described  by  Bard  (1974, 
p.  112-113).  Convergence  occurs  when  the  step  size  is  less  than  some 
specified  tolerance. 

To  implement  the  Newton-Greenstadt  algorithm  it  was  necessary  to 
evaluate  the  gradient  and  Hessian  for  c,  and  to  compute  the  eigenvalues 
and  eigenvectors  of  the  Hessian.  The  gradient  vector  was  computed 
directly  from  the  analytical  profile  expression;  the  Gauss  approximation 
to  the  Hessian  (Bard,  1974,  p.  96)  was  used.  This  latter  approximation 
was  chosen  solely  for  convenience  and  has  no  effect  on  the  accuracy  of 
the  final  result  since  the  Hessian  is  used  only  to  indicate  the  proper 
search  direction.  In  fact  at  any  stage  it  is  only  necessary  to  have  an 
approximation  to  H-.  The  eigenvalues  and  eigenvectors  of  H^  were  com- 
puted by  two  Purdue  University  Computer  Center  routines. 

In  practice  the  Newton-Greenstadt  procedure  performed  very  well. 

The  Hessian  H^  was  always  positive  definite.  Final  gradient  values  on 
the  order  of  10"®  and  smaller  were  typical.  However  in  some  cases, 
particularly  those  with  large  adverse  pressure  gradients,  it  was  neces- 
sary to  introduce  penalty  function  techniques.  Following  Bard  (1974, 
p.  159)  the  objective  function  was  modified  according  to: 

N 

G = e + I a./h.(x).  (2.2) 

i = l ’ ^ 

Here  e is  the  root-mean-square  error  as  previously  defined  and  h^(x) 
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represents  the  constraints  on  the  profile  parameters  such  that: 

h^.  (x)  > 0,  i = 1,  . . . , N . (2.3) 

Here  the  a/s  are  small  positive  constants  whose  values  are  dictated  by 
the  range  of  values  of  h^.  and  e (Bard,  1974).  The  gradient  vector  and 
Hessian  were  appropriately  altered  with  the  Gauss  approximation  applied 
to  H^ . The  basic  concept  behind  the  penalty  function  approach  is  to 
force  the  profile  parameters  away  from  values  which  they  cannot  assume; 
for  example  in  the  present  problem  u^  must  always  be  positive  and  if  a 
situation  was  encountered  in  the  optimization  procedure  where  u^  was 
approaching  zero,  then  the  procedure  given  in  equation  (2.2)  was  imple- 
mented with  say,  hj  = u^.  Once  the  minimum  of  the  redefined  objective 
function  in  equation  (2.2)  was  found,  a-^  was  reduced  by  a factor  of  10 
and  the  optimization  procedure  continued.  By  continuing  this  process 
iteratively  e is  eventually  reduced  arbitrarily  close  to  e in  equation 
(2.2).  This  technique  performed  very  well  in  the  more  difficult  cases. 

2.2  Direct  Search 

To  optimize  the  composite  profile  consisting  of  the  unsteady  wall 
layer  model  and  the  outer  similarity  profile,  a direct  search  minimiza- 
tion procedure  was  used.  This  technique  has  the  distinct  advantage  of 
only  requiring  values  of  the  objective  function.  This  is  important  in 
the  case  of  the  similarity  profile  since  gradients  of  the  objective 
function  (which  are  required  in  the  Newton-Greenstadt  procedure)  would 
have  to  be  evaluated  numerically.  The  basic  procedure  was  to  vary  one 
parameter  such  that  the  objective  function  always  decreased  while  hold- 
ing the  other  parameters  constant.  When  a local  minimum  had  been 
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bracketed,  quadratic  interpolation  was  used  to  refine  the  minimum 
location. 

More  specifically  an  initial  estimate  of  the  optimum  profile 
parameter  values  was  made  along  with  initial  parameter  step  sizes. 

The  objective  function  was  evaluated  at  the  initial  location.  One 
parameter  was  incremented  by  the  initially  specified  step  size  while 
the  remaining  parameters  were  held  fixed.  If  the  objective  function 
decreased,  then  the  stepping  process  was  repeated  until  the  objective 
function  increased;  this  bracketed  the  minimum.  If  the  objective 
function  initially  increased,  then  the  original  estimate  of  the  param- 
eter was  decreased  by  the  initially  specified  step  size.  As  before 
the  remaining  parameters  were  held  constant  and  the  process  was  con- 
tinued until  the  minimum  was  bracketed. 

Once  a minimum  had  been  bracketed,  a quadratic  interpolation  scheme 
was  applied  to  refine  the  parameter  value  at  the  minimum.  A quadratic 
curve  was  fitted  through  the  three  points  surrounding  the  minimum.  The 
resulting  expression  was  differentiated  to  determine  the  location  of  the 
curve's  minimum  and  the  objective  function  was  evaluated  at  this  refined 
parameter  value.  Since  the  objective  function  variation  was  not 
quadratic,  it  was  necessary  to  continue  the  interpolation  process  until 
there  was  no  appreciable  change  in  either  the  objective  function  or  the 
parameter  value. 

The  direct  search  minimization  procedure  described  above  functioned 
very  well.  Typically  less  than  100  objective  function  evaluations  were 
required  for  one  and  two  parameter  fits.  However  when  more  than  two 
parameters  are  optimized,  the  procedure  does  become  very  time  consuming. 
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2.3  Discussion 

For  both  the  Newton-Greenstadt  and  direct  search  minimization 
techniques,  considerable  effort  was  taken  to  assure  that  a valid  minimum 
had  been  found  for  each  velocity  profile  fit.  After  a candidate  best- 
fit  was  obtained,  different  initial  parameter  estimates  were  made  to 
determine  if  the  same  minimum  was  achieved.  Different  step  sizes  and 
error  criteria  were  also  used.  For  the  Newton-Greenstadt  procedure  the 
results  presented  in  this  report  are  for  a step  size  tolerance  of  10"®. 
However  essentially  the  same  fits  were  obtained  when  this  was  decreased 
to  10“^®.  For  the  direct  search  technique  the  relatively  small  step 
sizes  and  error  criteria  that  were  used  produced  only  very  slight  differ- 
ences in  the  fits.  In  either  case,  the  results  given  in  this  report  are 
those  with  the  smallest  root-mean-square  error,  c. 
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3.  INNER  REGION  DATA  COMPARISONS 

The  data  comparisons  between  the  Van  Driest  and  unsteady  wall 
layer  inner  region  mean  velocity  profile  models  will  be  presented  in 
this  section.  The  reason  for  considering  inner  region  comparisons  is 
two-fold.  First,  and  foremost,  they  provide  a direct  test  of  the  two 
inner  region  models  to  be  considered  without  introducing  any  question 
of  which  outer  region  profile  is  the  best.  A second  and  mainly  prac- 
tical reason  is  that  the  results  of  the  inner  region  optimization  serve 
as  good  initial  estimates  for  the  full  profile  fits  to  be  presented  in 
the  next  section. 

To  determine  those  points  which  truly  lie  within  the  wall  layer 
region,  the  experimentally  reported  u^  value  was  used  to  graph  lT''"  versus 
y'*'  on  the  standard  semi -log  plot.  The  inner  region  was  considered  to 
extend  from  the  wall  up  to  and  including  the  logarithmic  region;  however 
a practical  difficulty  exists  in  deciding  which  data  points  lie  in  the 
logarithmic  region  and  which  lie  in  the  outer  region.  The  exclusion  of 
data  points  which  are  considered  to  be  outer  region  points  is  somewhat 
subjective  but  in  practice  the  uncertainty  usually  concerns  one  and  at 
most  two  points. 

Before  moving  directly  to  a discussion  of  the  inner  region  data 
comparisons,  a description  of  the  two  inner  region  velocity  profile 
models  will  be  presented. 
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3.1  The  Van  Driest  and  Modified  Van  Driest  Profiles 

The  Van  Driest  model  was  originally  proposed  for  zero  pressure 
gradient  flows  (Van  Driest,  1956).  However,  because  of  its  apparent 
success  it  has  been  extensively  modified  to  account  for  non-zero  pres- 
sure gradients,  and  heat  and  mass  transfer.  Many  of  these  modifications 
are  discussed  by  Cebeci  and  Smith  (1974).  It  is  perhaps  unfair  to 
associate  the  name  Van  Driest  with  all  of  these  modifications  but  this 
is  common  practice  in  the  literature  and  thus  it  is  convenient  to  do  so 
in  this  report.  The  modified  model  forms  a basis  for  many  prediction 
methods  including  Cebeci-Smith,  Patankar  and  Spalding,  and  Bradshaw  and 
Ferris;  these  are  detailed  in  Kline,  et  al.  (1969). 

In  spite  of  the  many  variations,  the  basic  concept  of  the  Van 
Driest  model  remains  unaltered.  Prandtl's  mixing  length  is  used  to 
model  the  Reynolds  shear  stress  term  in  (1.2),  Specifically: 


uV  = 


8u^ 


3y* 


3y* 


Here  l is  a mixing  length  function  which  is  written  for  the  inner 
region  as: 


I - icy^F(y''’,T^,A'’') 


where 


F(y'*',T'‘’,A'^)  = 1 - exp 


m- 


(3.1) 


(3.2) 


Here  = t/t  is  the  total  shear  stress  and  A^  is  referred  to  as  the 

W 

viscous  damping  length  parameter.  The  total  shear  stress  can  be 
written  as 
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(3.3) 


where  u'*’  = u*/u^(x,Re).  If  equation  (3.3)  is  solved  for  3lTV9y^  and 
integrated  from  y^  = 0 to  an  arbitrary  point  within  the  wall  layer,  it 
follows  that 


+ Pt'*' 

U = ========  d^,  (3.4) 

K 1 + 

where  the  capital  notation  has  been  introduced  to  denote  the  portion 
of  u"*"  which  is  strictly  an  inner  region  profile.  It  follows  from 
integrating  the  momentum  equation  (1.2)  and  neglecting  the  convective 
terms  that  the  total  shear  stress  is 

= 1 + pV  (3.5) 


where 


V dUe 

— u — - 

jj3  e dx* 


(3.6) 


Equation  (3.5)  is  often  referred  to  as  the  Couette  approximation 
(Patankar  & Spalding,  1968)  and  it  is  of  interest  to  briefly  examine 
in  what  sense  (3.5)  is  an  approximation.  For  pipe  or  channel  flow 
equation  (3.5)  is  exact;  for  boundary-layer  flow  neglecting  the  con- 
vective terms  does  incur  an  error  in  (3.5)  for  any  finite  value  of  the 
Reynolds  number.  However  in  the  limit  Re  <»,  the  asymptotic  analysis  of 
Fendell  (1972)  suggests  that  to  leading  order  equation  (3.6)  is  correct. 

Substitution  of  (3.5)  and  (3.6)  into  (3.4)  results  in  an  expression 
for  U'*’  which  is  a function  of  the  two  parameters  A'*’  and  u^.  For  zero 
pressure  gradient  flow  Van  Driest  (1956)  found  that  a value  of  A'*'  = 26 
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represented  inner  region  data  reasonably  well;  when  (3.1)  is  used  in 
prediction  programs,  various  forms  are  used  for  A"*"  to  attempt  to  account 
for  pressure  gradient  effects  and  heat  and  mass  transfer.  Because  A'*' 
is  not  a physically  measureable  quantity  the  proposed  forms  for  A'*’  are 
often  obtained  indirectly  in  the  following  way.  A correlation  for  A'*' 
is  proposed  with  various  adjustable  constants  which  are  then  altered 
in  a prediction  program  so  that  profile  data  are  reproduced  as  well  as 
possible.  Of  course  there  are  many  sources  of  error  in  such  a procedure; 
the  alternative  that  is  adopted  here  is  that  both  u^  and  A'*’  are  regarded 
as  parameters  and  these  are  adjusted  to  best  fit  profile  data.  The 
result  of  this  procedure  is  the  best  representation  of  the  data  that  can 
be  obtained  by  equation  (3.4);  the  object  of  this  is  to  gain  some 
insight  into  how  well  the  Van  Driest  and  modified  Van  Driest  profile  can 
represent  the  inner  region  data. 

Because  the  integral  in  (3.5)  cannot  be  determined  analytically, 
for  a given  value  of  A'*’  and  u^  and  at  any  stage  in  the  optimization  pro- 
cedure, the  profile  was  calculated  numerically.  The  numerical  quadrature 
used  was  one  based  on  Simpson's  rule  and  to  achieve  a good  level  of 
accuracy  500  intervals  were  taken  between  the  wall  and  the  first  data 
point  and  then  250  intervals  between  subsequent  data  points.  Although 
this  integration  procedure  is  somewhat  time  consuming,  the  main  motiva- 
tion here  was  to  obtain  at  least  four  significant  figures  of  accuracy 
in  order  to  ensure  a fair  comparison. 

At  this  stage,  it  is  useful  to  point  out  that  in  spite  of  its  wide 
use  in  connection  with  non-zero  pressure  gradient  flows,  there  is  an 
inherent  contradiction  between  the  extended  Van  Driest  model  and 
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experimental  observation.  First  consider  the  case  of  zero  pressure 
gradient  flow  and  formally  take  p'*'  = 0 in  (3.5);  a simple  asymptotic 
expansion  of  the  derivative  of  (3.4),  for  large  y"*^  yields 

^ |l ^ - •••!  as  y'*’  ->■  00,  for  p'^  e 0.  (3.7) 

9y  <y  ^ Rk-  V 


2<y  8<  y 


Consequently  for  p^  = 0,  equation  (3.7)  implies  a logarithmic  behavior 
for  as  y^  becomes  large.  This  is  consistent  with  all  experimental 
measurements  of  the  time-mean  turbulent  velocity  profile  for  flows  at 
constant  pressure.  Thus  the  Van  Driest  (1956)  model  performs  correctly 
under  the  conditions  for  which  it  was  originally  proposed.  However  for 
any  non-zero  value  of  p"*"  a similar  analysis  yields; 

-1 


2p'^y‘''  2K^y‘*'»'p^^ 


as  y'*'  -*■«»,  p'*’  0. 


(3.8) 


Here  p"^  is  assumed  to  be  no  larger  than  0(1),  although  perhaps  small. 
Evidently  from  (3.8)  the  time-mean  profile  behaves  as  for  large  y'*’ 
and  there  is  no  logarithmic  term  in  (3.4)  either  to  leading  order  or  to 
higher  orders.  However  experiment  confirms  the  existence  of  a logarith- 
mic behavior  for  flows  with  non-zero  pressure  gradients.  This  fundamen- 
tal contradiction  with  well -documented  experimental  observation  of 
logarithmic  behavior  provides  a valid  reason  for  rejecting  the  modified 
Van  Driest  inner  region  model.  The  inability  of  the  modified  Van  Driest 
profile  to  properly  conform  to  flows  with  pressure  gradient  will  be 
borne  out  in  the  subsequent  data  comparisons.  First  however  a brief 
description  of  the  unsteady  wall  layer  model  will  be  given. 


3.2  Unsteady  Wall  Layer  Model 

The  unsteady  wall  layer  model  incorporates  the  observations  of 
several  inner  region  flow  visualization  studies  (Kline  & Rundstadle'r, 
1959;  Kline  et  al . , 1967;  and  Corino  & Brodkey,  1969)  with  an  analysis 
of  the  possible  motions  within  the  wall  layer  in  the  asymptotic  limit 
Re  Some  main  features  of  the  model  were  proposed  by  Black  (1968) 

based  partially  on  the  work  of  Einstein  and  Li  (1956).  Many  of  the 
assumptions  behind  Black's  (1968)  analysis  and  a number  of  features  of 
his  solution  for  the  wall  layer  flow  are  not  satisfactory.  These  are 
discussed  in  Walker,  Abbott  and  Scharnhorst  (1976)  and  Walker  and  Abbott 
(1977)  where  a complete  discussion  of  the  ideas  behind  the  unsteady  wall 
layer  model  is  given.  Only  a brief  discussion  will  be  given  here. 

It  is  observed  from  experiment  that  the  wall  region  flow  is  both 
three-dimensional  and  time-dependent  but  does  possess  a considerable 
degree  of  ordered  structure.  In  particular  two  features  dominate  the 
wall  layer  flow.  The  first  of  these  is  the  bursting  phenomenon  which 
can  be  regarded  as  a localized  breakdown  of  the  wall  layer  flow;  it  is 
characterized  by  rapid  and  violent  ejection  of  wall  layer  fluid  into 
the  upper  layer  and  is  of  relatively  short  duration.  This  event  is 
followed  almost  simultaneously  by  inrushing  fluid  from  the  outer  layer 
and  this  has  been  termed  the  sweep  by  Corino  and  Brodkey  (1969).  There 
then  ensues  a relatively  long  period  of  time,  called  the  quiescent 
period.  This  is  the  second  important  feature  of  the  wall  layer  flow 
and  it  is  during  this  period  of  time  that  the  streak  structure  is 
observed.  The  streaks  are  normally  visualized  with  hydrogen  bubble 
wire  techniques  and  the  tendency  of  the  bubbles  to  collect  into  streaks 
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indicates  surfaces  across  which  there  is  no  spanwise  motion.  Experi- 
ments indicate  that  the  average  spanwise  spacing  of  the  streaks, 

X"*”  = Xu^/v,  is  a number  on  the  order  of  100;  however  most  of  the 
experiments  that  suggest  a constant  value  of  X'*’  = 100  are  in  a fairly 
limited  Reynolds  number  range  and  the  high  Reynolds  number  results  of 
Gupta,  et  al.  (1971)  suggest  that  X'''  may  become  large  as  Re  «>.  In 
either  case.  Walker  and  Abbott  (1977)  have  argued  that  during  the 
quiescent  period  (when  the  integrity  of  the  wall  layer  is  intact)  the 
equations  governing  all  three  velocity  components  in  the  wall  layer 
reduce  to  linear  equations  of  the  heat  conduction  type.  The  consequence 
of  this  is  that  the  wall  layer  problem  during  the  quiescent  period  is  a 
linear  one  and  the  principle  of  superposition  of  solutions  applies. 
Between  two  adjacent  streaks  the  solutions  for  the  velocity  components 
may  be  written  as  a Fourier  series  in  z'*’  = z/X;  however  over  a long 
period  of  time  and  after  a considerable  number  of  cycles  only  that 
component  of  u"^  = u(x^,y’^,z'*',t''’)/u^  which  is  independent  of  z"*"  during 
each  quiescent  period  can  contribute  to  the  mean  profile.  Essentially 
this  is  because  after  each  burst  and  subsequent  quiescent  period  the 
streaks  are  observed  to  rearrange  in  different  z'*'  locations  and  over  a 
long  period  of  time  the  distribution  appears  random.  In  averaging  the 
time-dependent  solutions  for  u'*’  over  a single  period,  the  average  is  to 
be  interpreted  essentially  as  an  average  over  a large  number  of 
quiescent  periods.  Moreover  it  is  also  argued  (Walker  & Abbott,  1977) 
that  the  profile  so  obtained  represents  a good  approximation  to  the 
total  mean  profile  U'*'  in  the  wall  layer.  The  reason  for  this  is 
associated  with  experimental  observation;  if  attention  is  focused  on 
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a fixed  small  area  of  the  plate,  the  wall  layer  is  observed  to  be  in  the 

quiescent  state  in  excess  of  ninety  percent  of  the  total  time.  Thus  it 

may  be  that  in  the  limit  Re  -»■  <», the  time  scale  associated  with  the 

bursting  itself  is  small  with  respect  to  the  time  scale  associated  with 

the  quiescent  period.  The  main  difficulty  here  is  that  to  date  there 
is  no  argument  which  remotely  explains  the  presence  of  the  streaks  or 
the  burst.  The  analytical  difficulties  in  resolving  these  questions 
are  formidable  and  at  this  stage  neglecting  the  contribution  to  the 
mean  profile  from  the  burst  itself  is  conjectural  and  motivated  pri- 
marily from  experiment. 

In  Walker,  Abbott  and  Scharnhorst  (1976)  an  initial  condition  for 
u‘‘’(x''’,y'*',z'*',f^)  was  assumed  which  was  selected  to  model  the  latter 
stages  of  the  sweep.  The  following  expression  for  the  mean  profile 
was  obtained^: 


U'^  = -^  (1  + Yj  + 2-^ogHj)  I ^ He"^^ 


Here 


H 


+ HE*  (H)  + ^ cA|(H  + 2H2|e(H)  - ^ loqW  - ^ 


2p+$2  i 


+ |2h2£c-9(H/HJ  - ^ (H  - HJ(2H^  + 5HHi  - H^jj 

(3.9) 


= yV2S,  = y-;/2S,  S = 


t In  the  expressions  below  and  throughout  this  report  fogx  denotes  the 
natural  logarithm. 


r 


and  T is  the  mean  period  between  bursts;  L is  the  Heaviside  function 
and  Yq  is  Euler's  constant.  In  addition  the  function  H appearing  in 
(3.9)  is  defined  as 

-72  w2  (-y  _„2 

5(n)  = e ^ e-^  e ^ dxdydz.  (3.1 

■^0  "^0  "^0 


A partial  list  of  the  properties  of  (3.10)  is  given  in  Walker,  et  al. 
(1976)  and  a subsequent  report  in  ..his  series  will  give  a complete 
compilation. 

The  parameter  y|  appears  in  the  assumed  functional  form  of  the 
initial  condition;  it  can  be  related  to  the  log-law  constant  and  S 
from  two  profile  compatibility  relations  as  follows.  First,  upon 
expanding  (3.9)  for  large  y'*’  it  follows  that 

C^-  = ^ (1  - iog  yt)  - I p-^S^.  (3.11 

Second  from  the  definition  of  the  time-mean  velocity  profile 

— =1.  (3.1Z 

and  it  follows  from  (3.9)  and  (3.12)  that 

^ kS  + I p\s'  = -logH,  - ^ ^ ^ h!;  + 0(hJ). 


Equations  (3.11)  and  (3.13)  can  be  used  to  determine  values  of  , S 
and  y|  given  any  one  value;  in  practice  it  is  the  quantity  which  is 
usually  assigned.  Thus  equations  (3.9),  (3.11)  and  (3.13)  define  the 
unsteady  wall  layer  time-mean  velocity  profile  model  for  the  inner 
region  of  a two-dimensional  turbulent  boundary  layer. 


22 


In  Walker,  et  al . (1976)  two  objections  to  the  profile  (3.9)  were 
discussed;  first  the  profile  has  a jump  discontinuity  in  the  fourth 
derivative  at  y|.  This  was  shown  to  be  a direct  consequence  of  the 
assumed  initial  condition;  although  the  difficulty  is  of  little  conse- 
quence insofar  as  data  comparisons  are  concerned  (because  y’J'  is  in 
general  small),  it  is  clearly  unsatisfactory  on  a theoretical  basis. 
The  second  difficulty  is  that  (3.9)  does  not  satisfy  the  wall  compati- 
bility condition 


3||+ 


9"U 

3y+' 


= 0. 


y'^=0 


(3.14) 


This  condition  can  be  demonstrated  in  a variety  of  ways  (Rotta,  1962). 

A more  satisfactory  expression  for  U'*'  may  be  obtained  by  means  of  a 
similarity  analysis.  The  details  of  this  are  involved  and  are  given  by 
Walker  and  Scharnhorst  (1977);  only  the  final  result  will  be  quoted 
here.  The  main  idea  is  to  write  the  time-dependent  normal  coordinate 
as: 


n = 


2*'t+  + t! 


(3.15) 


rather  than  the  corresponding  coordinate  selected  in  Walker,  et  al . 
(1976).  The  motivation  for  the  introduction  of  t^  in  (3.15)  is  sug- 
gested by  the  rapid  character  of  the  wall  layer  breakdown  (the  burst) 
and  almost  immediate  re-establishment  of  the  streak  structure;  the 
non-linearities  that  are  almost  certainly  present  in  the  governing 
equations  during  the  burst  then  must  rapidly  subside  as  the  quiescent 
period  commences  at  t'*'  = 0 (by  definition).  In  effect  the  parameter  t„ 
in  (3.15)  corresponds  to  an  uncertainty  in  the  origin  of  time. 
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An  analysis  described  in  Walker  and  Scharnhorst  (1977)  leads  to 
the  result: 

= [l  + ||]|R(S^t^)Q(H)  + Z(H)  + P(S^t■J)W(H)] 

|R(0,t;)Q(H„)  + Z(H„)  + P(0,t+)W(Ho^, 


It 

$2 


(3.16) 


where 
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R(S^t+)  E A„  + ^ log(S^  + tt). 
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(3.18) 
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5l  -H= 


- 3H^ 


(3.19) 

(3.20) 

(3.21) 


It  can  be  shown  tfiat  for  large  y"*",  the  asymptotic  form  of  the  profile 
given  by  (3.16)  is 


~ i togy^  + C. 


(3.22) 


where  the  constant  A^  in  (3.17)  is  related  to  C^.  by 


(3.23) 


T 
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In  addition  it  is  insisted  that  (3.12)  and  (3.14)  be  satisfied  and 
\ this  leads  to: 


J 


($2  + to)‘^^[R(S2,t^)  + 3P(SStti] 

- (to)"^^[R(0,t^  + 3P(0,t^)]  = 0,  (3.25) 

respectively.  Consequently  with  (3.23)  the  three  parameters  C^. , S,  and 
to  in  (3.16),  can  be  determined  given  any  one  value,  using  (3.24)  and 
(3.25).  The  profile  (3.16)  may  easily  be  evaluated  numerically  and 
involves  approximately  the  same  amount  of  execution  time  as  is  necessary 
to  perform  the  numerical  integration  in  (3.4)  to  evaluate  the  Van  Driest 
profile.  A computer  program  to  evaluate  (3.16)  is  available  from  the 
authors  upon  request. 

The  profile  (3.16)  is  the  preferred  profile  since  the  objections 
discussed  by  Walker,  et  al.  (1976)  to  (3.9)  are  removed.  However  for 
all  practical  purposes  both  profiles  are  essentially  the  same.  In 
Figure  1,  both  profiles  are  plotted  for  C-j  = 5.0  and  k = 0.41  and  appear 
as  a single  curve.  For  (3.9),  equations  (3.11)  and  (3.13)  yield  S = 
10.52  and  y|  = 0.350  while  (3.24)  and  (3.25)  give  S = 10.50  and  t^  = 

8.01  X 10"\  Although  there  is  a slight  difference  in  S the  profiles 
themselves  are  identical  in  most  cases  to  three  significant  figures  as 
demonstrated  in  Table  1.  In  this  report  the  inner  region  fits  and  the 


Table  1 


Unsteady  Wall  Layer  Profile  Comparisons 


+ 

y 

Eq.  (3.9) 

Eg.  (3.16) 

0.5 

0.4999 

0.4992 

1 

0.9934 

0.9924 

2 

1.9494 

1.9478 

4 

3.7111 

3.7087 

6 

5.2609 

5.2579 

8 

6.6055 

6.6018 

10 

7.7607 

7.7567 

15 

9.9521 

9.9474 

20 

11.3965 

11.3914 

25 

12.3563 

12.3512 

30 

13.0199 

13.0146 

40 

13.8913 

13.8860 

50 

14.4858 

14.4805 

60 

14.9513 

14.9460 

80 

15.6714 

15.6661 

100 

16.2236 

16.2184 

t Because  of  the  algebraic  decay  in  both  (3.9)  and  (3.16)  the 
profiles  do  not  agree  precisely  at  y"^  = 100,  for  example. 
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full  profile  fits  with  Coles  "law  of  the  wake"  were  carried  out  with 
(3.9);  the  improved  inner  region  profile  (3.16)  was  used  with  the  outer 
region  similarity  profile  (Fendell,  1972)  for  the  full  profile  fits  to 
be  described  in  Section  5.  Since  it  has  been  demonstrated  in  Table  1 
that  for  practical  purposes  the  two  inner  region  profiles  (3.9)  and 
(3.16)  are  essentially  equivalent,  the  fact  that  two  different  forms 
of  the  unsteady  wall  layer  model  were  used  should  have  no  effect  on  the 
final  results  and  conclusions. 

Comparisons  with  inner  region  data  were  made  with  the  data  sets 
grouped  according  to  pressure  gradient.  The  three  non-transpiring 
flows  of  Andersen,  Kays  and  Moffatt  (1972)  were  considered;  these  flows 
are  at  zero,  mild  adverse  and  strong  adverse  pressure  gradients,  respec- 
tively, and  are  particularly  attractive  data  sets  because  there  are  a 
considerable  number  of  inner  region  data  points.  Comparisons  were  also 
made  with  the  Fraser  Flow  A from  the  Stanford  Conference  (Coles  & Hirst, 
1969).  Although  there  are  not  as  many  inner  region  data  points,  the 
pressure  gradient  is  more  severe  than  in  the  Andersen,  et  al.  (1972) 
flows. 

The  detailed  comparisons  are  presented  in  tabular  and  graphical 
form  in  Appendix  A,  to  which  the  reader  is  referred.  For  convenience, 
a detailed  discussion  of  the  results  for  each  data  set  is  also  included 
in  Appendix  A. 

3.3  Summary 

The  four  inner  region  data  comparisons  are  discussed  in  detail  in 
Appendix  A and  are  summarized  in  Table  2.  It  is  seen  from  the  values  of 
the  average  root-mean-square  error,  e,  that  the  extended  Van  Driest 


Flow 

Van  Driest 

Unsteady  Model 

e 

avg 

e 

S 

avg 

Zero  P.G. 
8100 

0.0072 

31.6 

0.0073 

10.70 

Mild  P.G. 
8200 

0.0076 

- 

0.0058 

- 

Strong  P.G. 
8300 

0.0070 

- 

0.0062 

- 

Fraser  A 
5000 

0.0181 

0.0066 

. 

model  fails  to  conform  to  velocity  profiles  in  adverse  pressure  gradi- 
ents. In  view  of  the  analysis  of  §3.1,  and  particularly  equation  (3.8), 
the  result  is  not  surprising.  Current  boundary-layer  prediction 
schemes  do  not  do  particularly  well  in  predicting  flows  where  the 
pressure  gradient  is  important;  one  reason  to  expect  difficulty  from 
predictive  schemes  that  use  a modified  Van  Driest  model  as  a basis  for 
the  wall  region  is  evident  from  (3.8).  It  is  thus  reasonable  to 
suggest  that  the  modified  Van  Driest  model  should  be  abandoned, 
although  it  is  so  firmly  entrenched  that  this  suggestion  is  not  likely 
to  be  readily  accepted.  It  should  also  be  noted  that  even  for  the 
zero  pressure  gradient  case  for  which  the  original  model  was  developed, 
the  arithmetic  average  value  varies  substantially  from  the  value 
A'*’  = 26.0  suggested  by  Van  Driest  (1956). 

Referring  again  to  Appendix  A and  Table  2 it  is  shown  that  the 
unsteady  wall  layer  model  exhibits  no  difficulties  in  fitting  inner 
region  data  in  flows  with  pressure  gradient.  In  fact,  it  is  seen  that 
with  increasing  adverse  pressure  gradient,  the  average  error  e remains 
essentially  unchanged  from  case  to  case.  It  should  also  be  noted  that  for 
the  zero  pressure  gradient  case,  the  arithmetic  average  value  is  in 

good  agreement  with  the  theoretical  value  of  10.50  discussed  in  §3.2. 

Finally  it  should  be  mentioned  that  the  four  data  comparisons 
presented  in  this  report  are  only  representative  of  those  obtained  with 
a large  number  of  other  flow  cases. 

In  the  next  two  sections,  the  unsteady  wall  layer  model  will  be 
combined  with  two  different  outer  region  models  to  form  a composite 
velocity  profile  over  the  entire  boundary  layer;  in  this  way  complete 
profile  data  may  be  compared; 


4.  FULL  PROFILE  DATA  COMPARISONS 
(LAW  OF  THE  WAKE) 

In  order  to  obtain  an  expression  for  the  velocity  profile  across 
the  entire  turbulent  boundary  layer  it  is  necessary  to  couple  an  inner 
region  profile  discussed  in  section  3 with  a model  for  the  outer  region. 
In  this  report  two  outer  region  models  will  be  considered;  in  this 
section  Coles'  (1956)  "law  of  the  wake"  will  be  considered  and  in  section 
5 a similarity  profile,  described  originally  by  Mellor  and  Gibson  (1966), 
will  be  investigated.  For  reasons  discussed  in  section  3,  the  modified 
Van  Driest  profile  for  the  inner  region  will  not  be  considered  further 
since  both  outer  region  profiles  considered  in  this  report  require  that 
the  inner  region  profile  behave  logarithmically  for  large  y'*'.  Also,  in 
view  of  the  results  in  section  3 both  outer  region  profiles  will  be 
coupled  with  the  unsteady  wall  layer  model  rather  than  the  Van  Driest 
(1956)  model. 

The  majority  of  the  incompressible  data  considered  here  is  taken 
directly  from  Coles  and  Hirst  (1969);  in  addition  the  more  recent  non- 
transpired  data  of  Andersen,  et  al.  (1972)  and  the  adverse  pressure 
gradient  data  of  Samuel  and  Joubert  (1974)  will  be  considered.  The 
values  of  mainstream  velocity  gradient  dUg/dx*  were  taken  directly  from 
Coles  and  Hirst  (1969)  and  from  the  appropriate  reports  for  the  more 


recent  data. 
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A major  objective  of  this  and  the  following  section  is  to  decide 
which  outer  region  model  best  represents  the  experimental  data;  in  addi- 
tion an  important  goal  is  to  decide  which  prpfile  model  is  best  suited 
for  a prediction  method.  These  points  will  be  considered  subsequently 
in  section  6. 


4.1  The  Law  of  the  Wake 

Coles'  (1956)  "law  of  the  wake"  is  a simple  formula  for  the  outer 
region  profile  and  is  an  empirical  correlation  based  on  the  observation 
of  data  in  a wide  variety  of  flows,  with  and  without  pressure  gradients. 
Coles  (1956)  postulated  that  the  velocity  deviation  in  the  outer  region 
from  the  law  of  the  wall,  when  normalized  with  respect  to  the  maximum 
deviation  at  the  boundary  layer  edge,  should  be  a single-valued  function 
of  y/6.  Here  6 is  the  boundary-layer  thickness.  When  several  data  sets 
were  reduced  to  this  form,  hey  appeared  to  exhibit  a common  behavior 
which  was  correlated  as  a numerical  wake  function  w(y/6).  It  was  demon- 
strated by  Hinze  (1959)  that  the  wake  function  could  be  closely  approxi- 
mated by 


w(y/6)  = 2 


(4.1) 


and  this  simple  analytical  expression  was  subsequently  used  by  Coles  and 
Hirst  (1969).  The  law  of  the  wake  was  then  proposed  as  a linear  combina- 
tion of  a logarithm  and  the  wake  function,  according  to. 


— * 

^0  1»  -2  firy' 

y +C.+— 


(4.2) 


Here  Uq  denotes  the  time-mean  velocity  profile  in  the  outer  region  of 
the  boundary  layer  and  II  is  the  so-called  wake  strength  parameter. 
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It  should  be  emphasized  that  (4.2)  is  not  a physical  law  in  the 
usual  sense  but  simply  a good  experimental  correlation.  The  main  focus 
in  this  section  will  be  on  how  well  (4.2)  represents  experimental  data 
when  coupled  with  the  unsteady  wall  layer  model;  we  return  to  a critical 
evaluation  of  this  model  in  section  6. 

If  K and  are  assumed  to  be  universal  constants,  equation  (4.2) 
contains  the  three  parameters  n,  6,  and  u^.  The  boundary-layer  thick- 
ness 6 is  normally  defined  as  that  physical  location  where  the  boundary- 
layer  profile  reaches  a specified  percentage  of  the  mainstream  value  with 
99%  and  99.5%  being  typical  experimental  definitions.  The  definition  of 
6 makes  this  quantity  a somewhat  nebulous  quantity  and  in  (4.2)  6 is 
generally  regarded  as  a parameter.  By  insisting  that  iJJ  = at  y = 6, 
where  is  the  local  mainstream  velocity,  the  so-called  skin  friction 
relation 


n = I ^ 


, 6u 

-a  - i tog  ^ 

U^  < ^ V 


- s- 


(4.3) 


is  obtained.  Equation  (4.3)  is  a relation  connecting  the  three  param- 
eters u^,  n and  6;  consequently  for  assumed  values  of  k and  , equation 
(4.2)  is  a two-parameter  expression  for  the  outer  region  profile. 

In  preparation  for  the  Stanford  Conference,  Coles  and  Hirst  (1969) 
assumed  that  the  values  of  k = 0.41  and  C.  = 5.0  were  universal  constants 
and  optimized  on  the  two  parameters  u^  and  6;  the  best  fit  to  the  data 
was  defined  as  in  equation  (1.5).  Both  values  of  k and  C^.  are  question- 
able and  other  values  are  often  quoted  in  the  literature.  The  main 
difficulty  in  deciding  whether  or  not  < and  C^  are  universal  constants 
is  that  there  is  as  yet,  no  direct  method  for  measuring  u^  in  a boundary- 


I 


1 


33 


layer  flow.  The  above  values  of  k and  are  essentially  extrapolated 
from  pipe  flow  data  where  can  be  measured  directly.  Nevertheless 
in  this  report,  a value  of  k = 0.41  will  be  assigned  tentatively. 

A composite  profile  consisting  of  the  unsteady  wall  layer  model 
and  the  law  of  the  wake  (4.2)  is  defined  by: 

u+=  + W(y/6),  (4.4) 


where 


W(y/6)  = 


2n 

K 


for  I < 1. 
for  y > 6. 


(4.5) 


The  additional  definition  in  the  second  part  of  (4.5)  is  made  so  that 
all  measured  data  points  can  be  reasonably  included  in  the  fitting  pro- 
cedure. In  Coles  and  Hirst  (1969),  data  points  near  the  mainstream  and 
in  some  cases  well  within  the  boundary  layer  were  omitted;  the  motivation 
for  this  was  that  the  slope  of  the  profile  defined  by  (4.2)  does  not 
approach  zero  as  y 6.  In  addition,  some  data  near  the  wall  were  also 
omitted  on  the  grounds  that  it  may  be  in  error.  The  rejection  of  data 
points  is  rather  subjective  and  in  this  report  all  measured  data  points 
were  included  in  the  comparisons. 

The  composite  profile  (4.4)  contains  the  three  independent  param- 
eters u^,  S and  ; one  main  difference  from  the  fitting  procedure 
carried  out  by  Coles  and  Hirst  (1969)  is  that  (4.4)  contains  an  adjust- 
able log-law  "constant."  At  any  stage  in  the  optimization  procedure, 
given  a value  of  C^. , values  of  t^  and  S may  be  determined  from  (3.23), 


i 
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(3.24)  and  (3.25)  and  the  current  form  of  is  completely  determined 
from  (3.16). 

The  data  fits  are  given  in  graphical  and  tabular  form  in  Appendix  B 
to  which  the  reader  is  referred.  The  data  sets  are  arranged  and  numbered 
in  the  same  order  as  in  Coles  and  Hirst  (1969).  The  last  two  data  sets 
of  Samuel  and  Joubert  (1974)  and  Andersen,  et  al.  (1972)  have  been 
labelled  as  7000  and  8000  series  respectively.  For  convenience,  brief 
comments  about  each  data  set  and  the  corresponding  fit  are  included  in 
Appendix  B. 
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5.  OUTER  REGION  SIMILARITY  PROFILE 


5.1  Introduction 

The  existence  of  self-similar  profiles  for  two-dimensional  time- 
mean  turbulent  boundary  layers  is  a subject  which  has  attracted  a good 
deal  of  attention  in  the  literature.  It  is  important  here  to  make  a 
distinction  between  the  terms  equilibrium  and  self-similar  which  are 
often  used  interchangeably.  The  term  equilibrium  is  understood  to  apply 
to  a turbulent  boundary  layer  for  which  the  mainstream  velocity  is  such 
that  either  U (x)  ~ x“  or  U (x)  ~ e^^  where  a and  B are  constants.  In 
order  for  self-similar  profiles  to  exist  equilibrium  is  a necessary  but 
by  no  means  a sufficient  condition. 

To  investigate  this  point  further,  consider  the  case  of  laminar 
boundary  layers  where  similar  solutions  satisfy  the  well  known  Falkner- 
Skan  equation.  The  existence  of  such  solutions  in  laminar  flows  has 
been  discussed  by  Brown  and  Stewartson  (1965)  who  argue  that  similar 
solutions  may  be  expected  in  two  physical  situations.  The  first  of  these 
is  at  an  x station  where  the  velocity  profile  is  an  initiator  of  the 
boundary-layer  flow  downstream;  such  a situation  occurs  physically  at 
the  front  stagnation  point  of  a bluff  body  or  the  leading  edge  of  a flat 
plate,  for  example.  In  these  cases  the  Falkner-Skan  profile  gives  the 
proper  laminar  boundary-layer  solution  at  a point  of  attachment  of  the 
mainstream  and  provides  the  initial  condition  to  initiate  a boundary- 
layer  calculation  downstream.  For  turbulent  boundary  layers  there  would 
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appear  to  be  no  analogue  of  this  physical  situation.  Generally  the 
flow  at  a point  of  attachment  of  the  mainstream  is  observed  to  be 
laminar  and  when  the  downstream  boundary-layer  flow  is  turbulent,  there 
is  a transition  zone  in  between.  Moreover,  a wide  variety  of  experi- 
mental conditions  can  lead  to  transition  and  an  eventual  fully 
developed  turbulent  boundary  layer  downstream  for  the  same  mainstream 
velocity  distribution.  Consequently,  there  would  appear  to  be  no 
reason  to  expect  that  the  flow  in  the  outer  region  of  an  equilibrium 
boundary  layer  will  be  self-similar  at  the  initial  stations  of  a fully 
developed  turbulent  flow. 

For  laminar  boundary  layers,  the  other  case  discussed  by  Brown  and 
Stewartson  (1965)  is  that  where  a similarity  solution  becomes  what  may 
be  described  as  a terminator  of  a more  general  boundary-layer  flow  and 
two  cases  of  this  behavior  are  considered  by  Brown  and  Stewartson  (1965). 

A physical  case  where  this  can  occur  is  at  a point  of  detachment  of  the 
inviscid  flow;  examples  of  such  behavior  are  known  in  magnetohydrodynamic 
flow  (Leibovich,  1967a,  1967b)  and  in  rotating  flows  (Walker  & Stewart- 
son, 1972)  at  the  rear  stagnation  point  of  symmetrical  and  two-dimensional 
bluff  bodies.  Another  case  is  flat  plate  flow  and  here  if  the  initial 
velocity  profile  at  any  arbitrary  location  on  the  plate  is  not  given  by 
the  Blasius  solution,  then  the  Blasius  profile  can  only  become  the 
relevant  solution  at  an  infinite  distance  downstream.  For  turbulent 
boundary-layer  flows  an  analogous  type  of  situation  is  expected;  that  is 
similarity  solutions  are  only  anticipated  as  termindtors  and  usually  at 
an  infinite  distance  downstream  from  wherever  the  boundary  layer  is 


created. 
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In  practice,  one  would  expect  to  be  able  to  measure  turbulent 
boundary-layer  profiles,  at  large  distances  downstream  of  the  transi- 
tion zone,  which  become  arbitrarily  close  to  being  self-similar;  how- 
ever, near  the  transition  zone  there  is  no  reason  to  expect  the  same 
behavior.  In  zero  and  favorable  pressure  gradients,  measured  profiles 
in  an  equilibrium  flow  should  increasingly  approach  self-similarity  at 
subsequent  data  stations  downstream.  In  adverse  pressure  gradient 
flows,  boundary-layer  separation  will  often  occur  before  similarity  is 
achieved. 


5.2  Similarity  in  the  Outer  Region 

An  ordinary  differential  equation  for  the  outer  region  similarity 
profile  has  been  obtained  by  a number  of  authors,  most  notably  Mellor 
and  Gibson  (1966),  Yajnik  (1970)  and  Fendell  (1972).  The  approaches 
used  by  these  authors  differ  somewhat  but  the  end  result  is  essentially 
the  same.  In  this  section,  the  spirit  of  the  approach  of  Fendell  (1972) 
will  be  adopted  (although  the  notation  is  slightly  different). 

Fendell  (1972)  considered  the  asymptotic  structure  of  the  time- 
mean  turbulent  boundary-layer  equations  in  the  limit  as  Re  ->-  « and 
showed,  on  the  basis  of  a few  reasonable  assumptions,  that  the  boundary 
layer  is  a composite  double  layer.  It  was  assumed  that  the  turbulence 
could  be  modelled  as  an  additive  stress  term  a so  that  the  time-mean 
momentum  equation  is. 


7 = u 

9y  ” dx 


9g 

9y 


1 

Re 


3^u 


(5.1) 


Here  all  velocities  and  lengths  have  been  made  dimensionless  with 
respect  to  Uq  and  L,  a representative  velocity  and  length  respectively 
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and  Re  = U^L/v  is  the  Reynolds  number.  For  practical  purposes  the 
function  a may  be  related  to  the  u'v'  product  by 


a = 


u'v' 


(5.2) 


Implicit  in  (5.2)  is  the  notion  that  the  intensities  u'^  and  v'^  and 
the  u'v'  product  are  of  the  same  order  of  magnitude;  consequently  the 
intensities  that  are  often  included  in  equation  (5.1)  may  reasonably 
be  neglected  since  they  are  differentiated  on  x and  will  be  of  lower 
relative  order. 

Fendell  (1972)  assumed  that  in  the  outer  region  the  velocity  pro- 
file may  be  written  as  a linearization  about  the  mainstream  value;  it 
can  be  shown  that  this  assumption  gives  rise  to  a consistent  asymptotic 
structure  and  in  fact  becomes  the  experimentally  determined  defect  law. 
For  the  outer  region, 

u*(x;Re)  (n,x)  + •••!  (5.3) 

o = U_^(x)|u*^(x;Re)  Z^{r\,x)  + •••!  (5.4) 

Here 

n = y/Ag(x;Re)  , (5.5) 


is  the  scaled  normal  coordinate  in  the  outer  region  and  is  the 
dimensionless  outer  region  length  scale.  In  addition,  u*  = u^(x;Re)/Ug(x) 
is  an  important  parameter  which  emerges  from  Fendell's  (1972)  analysis 
and  can  be  shown  to  be  0(1/ fog  Re^)  where  Re^  is  the  Reynolds  number 
based  on  the  local  boundary-layer  thickness;  also  it  emerges  that  Ag 
is  0(u*)  as  Re  -►  it  is  worthwhile  to  note  that  (5.4)  is  equivalent 
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to  assuming  that  u'v’  is  O(u^).  This  assumption  as  well  as  equations 
(5.2)  and  (5.3)  are  not  expected  to  apply  immediately  upstream  of  a 
point  of  time-mean  separation. 

The  equation  for  3Fj/3ri  is  obtained  by  substitution  of  (5.3)  and 
(5.4)  into  (5.1).  To  leading  order  the  laminar  shear  term  in  (5.1)  is 
negligible  and  the  dominant  terms  are  the  convection  and  Reynolds 
stress  terms.  Consequently 


00 


, Ao  3^Fi 
^ 3n^  u*'3n3x 


(5.6) 


with 


Fi(x,ri)  -»■  0 as  n 0, 
and 


(5.7) 


1^  (x,n)  -»■  0 as  q ->  <».  (5.8) 

In  equation  (5.6)  a prime  denotes  differentiation  with  respect  to  x. 

The  necessary  conditions  for  similarity  follow  from  (5.6)  and  are 

that 


= -26  = constant,  (5.9) 

^00 

and 

= a = constant.  (5.10) 


In  order  to  make  further  progress,  it  is  necessary  to  discuss  the 
nature  of  the  inner  region. 
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The  expansion  given  in  (5.3)  is  not  uniformly  valid  and  in  order 
for  the  no  slip  condition  to  be  satisfied,  an  inner  layer  is  required. 
In  the  inner  region  the  leading  terms  in  the  expansions  for  velocity 
profile  and  stress  are  (Fendell,  1972), 

lT  = u*U  (x)  — L (x,y^)  + •••  , (5.11) 

3y 

and 

a = u*^U^(x)  Oj(x,y'*')  + •••  , (5.12) 

where 

u L 

y+  = — y = Re  u*Ujx)y  , (5.13) 

is  the  conventionally  defined  inner  region  variable.  Matching  of  (5.3) 
and  (5.11)  takes  place  in  the  limits  y'*'  -*-  « and  n -*■  0 for  fixed  x;  in 
order  that  both  the  vorticity  and  velocity  profile  merge  smoothly,  both 
profiles  must  behave  logarithmically  according  to, 

1^  ~ ^ fog  n + Co(x)  , as  n^O,  (5.14) 

and 

y""  + C.(x)  , as  y^-co.  (5.15) 

Here  Cg  and  C^.  are  in  general  functions  of  x;  k is  analogous  to  the 
von  Kerman  "constant."  It  is  often  assumed  on  the  basis  of  extrapola- 
tion from  pipe  flow  data  that  k is  constant  and  equal  to  about  0.41. 
Whether  or  not  this  extrapolation  is  valid  or  not  is  conjectural;  from 
a mathematical  standpoint,  whether  or  not  < depends  on  x depends  on  the 


particular  constitutive  model  adopted  for  u'V.  One  additional  feature 
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worth  noting  is  that  it  also  follc-'S  from  the  analysis  (Fendell,  1972  ) 


Ej  -»•  1 as  n 0 « 


(5.16) 


Oj  -»■  1 as  y'*’ 


(5.17) 


where  again  flows  in  the  immediate  vicinity  of  a point  of  time-mean 


separation  are  excluded. 


As  yet  the  function  A.(x,Re)  has  not  been  specified  other  than  to 


remark  that  it  is  0(u*)  as  Re  ->  oo.  A relation  connecting  u*  and  h. 


follows  from  matching  (5.3)  and  (5.11),  using  (5.14)  and  (5.15)  and  is 


^co(x)}  + C.  - Co  . 


(5.18) 


To  leading  order  as  Re  «,  equation  (5.18)  implies  that 


log  Re.  ’ 


(5.19) 


which  is  a result  referred  to  earlier. 


We  are  now  in  a position  to  return  to  the  question  of  self-similarity; 


in  order  to  assess  the  magnitude  of  the  ratio  u*'/u*  which  appears  in 


(5.9)  and  (5.10),  equation  (5.18)  is  differentiated  with  respect  to  x. 


to  obtain. 


= (5.20) 


Using  (5.18)  to  eliminate  the  first  term  on  the  right  side  of  (5.20), 


it  follows  that. 


I 
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u*'  - <'  yi  j.  u*'  + !^+K.'(r  _r  ) + iefC  - C vl  (521) 

u^'ir  k\ao^u*  ^ 1-  SJ  + ^o>  j- 

Although  < may  in  general  be  a function  of  x,  it  must  become  constant 
for  a self-similar  flow  and  it  follows  from  (5.21)  that 


nJ^=  o(u*)  . 


(5.22) 


Thus  terms  containing  this  ratio  may  be  neglected  to  leading  order  in 
(5.9)  and  (5.10)  which  become,  respectively. 


q if  = -6. 


(5.23) 


q'  + ^ q = a. 


(5.24) 


Here  the  quantity  q is  defined  by  q = Aj^/u*.  Combining  these  two  equa- 
tions gives 


q'  = a + 8. 


(5.25) 


Consequently  only  two  types  of  mainstream  velocity  distributions  can 
lead  to'  self-similar  solutions  and  these  are; 


(a)  U (x)  = Di(x  - Xg)  * , for  a + 8 0 


(5.26) 


where 


(5.27) 


with  Dj  and  x^  arbitrary  constants. 


(b)  U„(x)  = D,e  , 


for  a + 8 = 0 
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where 


“2  = X 


S 


q = ^ = 
^ u* 


(5.29) 


with  Dj  and  Xq  arbitrary  constants. 

Equilibrium  boundary  layers  have  received  a good  deal  of  experi- 
mental attention  and  have  been  defined  in  this  context  as  flows  in 
which  the  dimensionless  velocity  defect  (Ug-u*)/u^,  when  expressed  as 
a function  of  y/6,  becomes  close  to  being  invariant  with  downstream 
distance.  This  type  of  behavior  is  expectec  mathematically  as  a solu- 
tion approaches  self-similarity.  Clauser  (1954,  1956)  suggested  that 
the  criterion  for  equilibrium  should  be  a constant  value  of  the  param- 
eter 3j,  defined  as 


6*  '^^e 

^c-fr^edlT* 

T 

Here  6*  is  the  dimensionless  displacement  thickness  defined  as 


(5.30) 


6*  = 


(5.31) 


In  view  of  (5.3)  and  (5.5),  to  leading  order. 


6*  = - Ag  u*  Fj(x,oo)  + 


(5.32) 


As  a condition  of  self-similarity  is  approached  Fj(x,<»)  must  approach 
a constant  value,  say  Fj(«>);  thus  equation  (5.30)  becomes 


g = A^Li.(°°.).  ^ = _bf  (oo) 


(5.33) 


t Because  (5,3)  is  a defect  law,  9Fj/9n  and  Fj(x,»)  are  in  general 
negative. 
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Clearly  both  types  of  mainstream  velocity  distributions  (5.26)  and 
(5.28)  give  a constant  value  of  which  again  is  a necessary  condition 
for  self-similarity. 

5.3  The  Eddy  Viscosity 

In  order  to  define  an  outer  region  profile  it  is  necessary  to 
select  a constitutive  relation  for  the  u'v'  product  in  the  outer  region. 
It  has  traditionally  been  assumed  that  there  exists  some  unknown  func- 
tional relation  between  the  u'v'  product  and  the  mean  velocity  field. 
This  view  will  also  be  adopted  here;  however,  it  must  be  recognized 
that  no  such  general  functional  relationship  has  ever  been  shown  to 
exist  and  for  the  present  any  such  hypothesis  must  be  regarded  as  no 
more  than  a useful  and  convenient  approximation.  The  major  difficulty 
in  resolving  this  point  is  that  the  time-dependent  dynamics  of  the 
outer  layer  are  not  understood.  From  experiment  it  has  become  increas- 
ingly evident  in  recent  times  that  the  burst-sweep  sequence  provides 
the  major  contribution  to  the  u'v'  product  for  both  the  inner  and  outer 
layer.  Consequently,  it  would  appear  necessary  to  analyze  a typical 
burst-sweep  sequence  in  order  to  either  demonstrate  or  deny  the  exist- 
ence of  a general  functional  relationship  between  u'v'  and  u.  However, 
the  burst-sweep  sequence  is  a time-dependent  viscous-inviscid  inter- 
action between  the  inner  and  outer  layer  caused  by  some  unknown  physical 
mechanism  in  the  outer  layer.  Moreover  even  if  the  cause  of  the  burst 
were  known,  the  mathematical  difficulties  associated  with  such  a problem 
are  formidable.  Consequently  for  the  present,  it  will  be  assumed  that 
some  constitutive  relation  does  exist. 
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The  particular  relation  that  will  be  adopted  here  is  that  of  an 
eddy  viscosity  model  for  the  outer  layer.  The  main  reasons  for  this 
choice  are  simplicity  and  the  relative  degree  of  success  that  this  type 
of  assumption  has  enjoyed  in  prediction  methods.  Indeed,  at  present  it 
is  far  from  clear  that  the  increasing  variety  of  much  more  involved 
second-order  closure  schemes  offer  any  significant  advantage  in  connec- 
tion with  predicting  turbulent  boundary-layer  flows.  A good  discussion 
of  the  eddy  viscosity  assumption  is  given  by  Saffman  (1970)  which  is 
quoted  directly  below: 

The  assumption  of  an  eddy  viscosity  has  been  severely 
criticized  by  some  writers,  who  also  regard  it  as 
harmful.  The  author  regards  this  view  as  extreme  and 
unjustified.  The  eddy  viscosity  hypothesis  is  equiva- 
lent to  retaining  the  leading  term  in  an  expansion  of 
the  (unknown)  functional  relation  between  Reynolds 
stress  and  the  mean  velocity  distribution. 

Consequently  for  the  outer  region,  define 


u'V 


a- 


From  equations  (5.4)  and  (5.5), 
■ A„UoU^(x)Lu*  • 


(5.34) 


(5.35) 


In  selecting  a functional  form  for  Eq,  the  same  type  of  model  that  is 
used  in  the  Cebeci-Smith  (1969)  and  Mellor  and  Herring  (1969)  prediction 
methods  will  be  adopted.  This  particular  type  of  model  is  discussed  by 
Mellor  and  Gibson  (1966)  and  is  credited  to  Clauser  (1954,  1956)  by 
these  authors. 

, One  main  feature  of  this  model  is  that  for  fixed  x and  in  the  outer- 

most regions  of  the  outer  layer,  Cq  is  selected  to  approach  a constant 
j value.  Specifically, 


\ 
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€q  ->  KUg(x)6*L. 


(5.36) 


Here  6*  is  the  dimensionless  displacement  thickness  defined  in  equation 
(5.31)  and  K is  an  empirically  determined  constant;  Cebeci  and  Smith 
(1969)  take  K equal  to  0.0168  while  Mellor  and  Herring  (1969)  choose 
0.016.  For  small  n,  equation  (5.36)  is  clearly  not  appropriate  in 
view  of  (5.16)  since  Zj  must  approach  one  as  n 0.  Thus  the  stress 
function  is  defined  according  to 


6*  , ^ 3"F, 

^ 3^^' 


(5.37) 


where  the  function  e(n)  must  be  such  that, 


e -►  K as  n 


(5.38) 


e -»■  KTi 


as  n 0. 


(5.39) 


From  (5.32),  to  leading  order 


= -Fi(x,«>) 


(5.40) 


and  therefore  (5.37)  and  (5.39)  become 


= -Fi(x,oo)  e(n) 


(5.41) 


-Fi(x,<»)e  ->■  Kn  as 


(5.42) 


In  general  the  functional  form  of  e(n)  is  not  known;  however  the  limit 
(5.38)  is  assumed  known  from  experiment  while  the  limit  (5.42)  is  known 
from  theoretical  considerations.  A simple  form  is  to  assume  two  straight 
line  variations  according  to. 


K for  n > rij , 

E = = ■ (5.4: 

-F,U,.)  f'”' 

where  Hi  = -KF j (x ,“)/<.  The  above  eddy  viscosity  formulation  is  that 
used  by  Mellor  and  Gibson  (1966)  and  Mellor  and  Herring  (1969)  who  in 
addition  select, 


6* 

An  = A = . 

0 c u* 


(5.44) 


The  reason  for  this  particular  choice  of  the  outer  region  length  scale 
is  motivated  by  (5.42)  which  assumes  a particularly  simple  form  in  this 
case;  in  particular,  Fj(x,«>)  = -1  under  this  selection.  In  Figure  2, 
the  eddy  viscosity  function  given  by  (5.43)  is  plotted  for  < = 0.41  and 
K = 0.0168;  for  comparative  purposes  it  will  be  convenient  to  use  the 
outer  region  length  scale  given  by  (5.44)  in  Figure  2. 

One  clearly  objectionable  feature  about  Cn;|(n)  is  that  the  first 
derivative  is  discontinuous  at  n = hi-  This  behavior  is  rather 
unpalatable  since  if  an  eddy  viscosity  function  does  exist,  it  should 
be  expected  to  behave  smoothly  and  continuously  for  all  n-  In  selecting 
a functional  form  for  e(n)  it  should  also  be  expected  that  e(n)  will 
monotonical ly  increase  from  the  linear  behavior  in  (5.42)  near  n = 0 to 
the  constant  value  in  (5.38)  for  large  n-  There  is  an  endless  variety 
of  functions  which  will  satisfy  the  above  requirements  and  although  a 
number  of  functions  were  tried,  the  results  for  only  two  will  be 
discussed  here. 
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dn  j 


n 


dUi 

dn 


= 0 


(5.49) 


where  B is  obtained  from  (5.33)  and  B = -6*U'/(u*^U  ).  Equation 

V t qq/  \ qq/  I 

(5.47)  may  be  integrated  across  the  outer  layer  to  obtain 


^ = 1 + 3B^ 
u*2  c 


and  consequently  (5.49)  becomes 


dn 


dUi 

dn  . 


+ 2B  Uj  = 0. 


(5.50) 


Equation  (5.50)  contains  three  parameters,  namely  u*,  K and  k.  Setting 
K = 0.0168  and  k = 0.41  completely  specifies  the  eddy  viscosity  and  for 
a given  value  of  u*,  the. solution  of  (5.50)  may  be  calculated. 

A composite  profile  for  the  entire  boundary  layer  consisting  of 
the  unsteady  wall  layer  model  and  the  solution  of  (5.50)  may  be  defined 
according  to. 


UqU  _ u 
u u*L)  (x) 

J oo'  ' 


+ Uj (n)  - ^ fog  n - Cg , 


(5.51) 


where  U'*'  is  defined  by  (3.16).  The  match  condition  (5.18)  is 

^ = ^tog  |6*ReUjx)}  + C.  - C„,  (5.52) 

and  is  related  to  the  inner  region  parameter  S through  (3.23).  At 
the  data  station  (x  = 0)  where  the  fit  is  to  be  performed,  take  the 
representative  velocity  U,  equal  to  the  local  mainstream  value  so  that 
U^(0)  = 1;  in  addition  at  the  data  station  being  fit  the  measured  dis- 
placement thickness  6*L  is  taken  from  experiment.  Thus  with  k,  K 
specified  and  the  displacement  thickness  taken  from  experiment. 
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equation  (5.51)  is  one  parameter  family  of  profiles  with  the  single 
parameter  being  u^. 

Because  Uj  behaves  logarithmically  for  small  n,  the  accurate 
numerical  solution  of  (5.50)  is  non-trivial  and  it  is  worthwhile  to 
briefly  discuss  how  this  was  carried  out.  Near  n = 0 the  solution  for 
Uj  was  written  as  a series, 

Uj  = i n + C()  + ))  n'^(a^  log  n + c^)  . (5.53) 

1 

Substitution  into  (5.50)  yields  a set  of  algebraic  equations  for  the 
a^  and  c^;  however,  because  the  boundary  condition  (5.8)  has  not  been 
applied  at  this  stage,  the  series  solution  (5.53)  is  not  uniquely 
determined  and  contains  the  arbitrary  constant  C^.  A value  was  first 
guessed  for  Cq  and  the  series  evaluated  at  a non-zero  value  of  p = 0^, 
say.  In  view  of  the  singular  behavior  of  (5.53)  near  n = 0,  the  value 
of  pg  must  be  selected  large  enough  so  that  finite  difference  approxi- 
mations to  (5.50)  may  be  made  with  good  accuracy  for  p > Oq;  on  the 
other  hand  p^^  must  be  small  enough  so  that  the  series  (5.53)  converges. 
In  practice,  a value  of  p^^  = 0.02  was  found  to  be  adequate  for  this 
purpose.  A solution  for  n i Oq  was  then  calculated  numerically  using 
a method  similar  to  that  described  by  Keller  (1969).  This  method  is  a 
boundary  value  technique  which  was  implemented  by  defining  a non- 
uniform  numerical  mesh  in  the  range  Pg  <.  p < L.  Here  Lira  large  but 
finite  value  of  p;  the  value  of  L,  on  which  condition  (5.8)  was  imposed 
as  an  approximation,  must  be  increased  until  there  is  no  significant 
change  in  the  computed  solution.  The  numerical  mesh  that  was  used 
employed  a telescoping  spacing,  wherein  the  mesh  spacing  was  selected 
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as  initially  small  and  uniform  near  n = tIq;  the  mesh  spacing  was  then 
successively  increased  by  a factor  of  1.02.  It  can  be  shown  that  the 
result  of  this  procedure  is  a set  of  tridiagonal  difference  equations 
which  are  easily  solved  by  a direct  method. 

Having  computed  a solution  with  one  arbitrarily  selected  value  of 
Cq,  another  solution  was  calculated  with  another  arbitrarily  selected 
value  of  C;,.  Both  solutions  computed  in  this  way  do  not  have  a continu- 
ous first  derivative  at  Po  but  may  be  combined  in  such  a way  that  both 
Uj  and  dUj/dri  are  continuous  at  p = Pg.  The  essential  features  of  this 
type  of  procedure  are  discussed  by  Walker  and  Stewartson  (1974). 
Approximately  350  mesh  points  v/ere  used  on  the  interval  p^^  < p < L and 
this  was  found  to  be  sufficient  to  guarantee  at  least  three  significant 
figures  of  accuracy  in  the  computed  numerical  solution. 

Before  moving  to  the  data  comparisons,  it  is  necessary  to  discuss 
one  further  point.  In  the  analysis  leading  to  (5.47)  it  has  been 
implicitly  assumed  that  the  inner  region  pressure  gradient  parameter 
p'*'  (defined  by  equation  (3.7))  is  small  as  Re  ->-<».  This  assumption  was 
not  necessary  in  section  3 and  the  inner  region  profile  given  by  (3.16)  is 
also  valid  when  p"^  is  0(1).  However,  for  the  outer  layer,  when  p'*'  is 
0(1)  the  structure  of  the  asymptotic  expansions  must  be  altered  and 
(5.3)  and  (5.4)  are  no  longer  appropriate.  It  is  worthwhile  to  briefly 
consider  the  physical  situation  when 


+ - ^ - 
~ L u ^ (x  ,Re)  dx 

T 

is  0(1)  and  this  occurs  for 
U7U„  = 0(Reu*^). 


1 X 

Re  u*^  U dx 


(5.54) 


(5.55) 


With  X fixed,  the  right  side  of  (5.55)  becomes  infinite  as  Re  and 
it  is  clearly  unacceptable  for  p"*"  to  become  0(1)  in  this  limit  process. 
However,  with  Re  large  and  fixed,  a situation  may  be  envisaged  where 
p'*'  becomes  0(1)  as  x approaches  a point  of  time-mean  separation, 
specifically  when  u = 0(Re  ).  When  this  situation  occurs  it 
appears  that  the  expansions  in  the  outer  layer  for  both  the  velocity 
and  additive  stress  in  (5.3)  and  (5.4)  respectively,  are  no  longer 
complete.  In  addition,  the  constitutive  relation  defined  by  (5.35) 
and  (5.45)  is  no  longer  adequate.  This  difficulty  has  been  briefly 
mentioned  by  Yajnik  (1970)  and  Fendell  (1972)  but  the  proper  way  to 
handle  the  near-separation  problem  is  not  resolved  at  present.  Conse- 
quently, in  this  and  the  following  section,  p'*’  is  regarded  as  a lower 
order  term  for  the  inner  layer  and  the  inner  region  profile  (3.16) 
was  taken  with  p'*'  formally  set  equal  to  zero.  (Note  that  the  log-law 
"constant"  in  (3.16)  is  still  variable  though  through  the  parameter 
S.)  In  view  of  the  above  discussion  and  that  in  §5.1,  attention  will 
be  restricted  here  to  flows  which  are  not  close  to  separation. 

The  data  comparisons  for  a number  of  equilibrium  flows  using  the 
composite  profile  defined  by  (5.51)  are  given  in  Appendix  C to  which 
the  reader  is  referred.  Detailed  descriptions  of  each  case  are  also 
given  in  Appendix  C and  only  the  main  results  will  be  summarized  here. 

Generally  the  one-parameter  (u.^)  fits  to  the  constant  pressure  and 
favorable  pressure  gradient  equilibrium  flows  give  good  representations 
of  the  experimental  data;  in  accordance  with  the  discussion  in  §5.1, 
the  representation  of  the  data  improves  at  the  downstream  stations  as 
the  flows  presumably  approach  self-similar  conditions.  In  the  adverse 
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equilibrium  flow  considered  in  Appendix  C,  the  representation  of  the 
data  is  fair  but  the  trend  to  self-similar  behavior  was  not  observed; 
this  is  the  expected  result  since  self-similarity  is  not  anticipated 
in  an  adverse  pressure  gradient  flow. 

5.5  An  Alternative  Eddy  Viscosity 

One  property  of  the  smoothly  varying  eddy  viscosity  formula  used 
in  §5.4  is  that  £5(0)  approaches  the  required  linear  behavior  (equation 
(5.42))  algebraically  as  n -»•  0.  An  alternative  formula  may  be  defined 
such  that  e(n)  approaches  the  linear  behavior  near  n = 0 exponentially 
qU'Ckly.  It  was  found  that  such  an  eddy  viscosity  relation  produced 
slightly  better  data  fits  and  consequently  this  form  of  the  eddy  vis- 
cosity will  be  briefly  considered  here. 

One  possible  choice,  satisfying  the  above  requirement,  would  be 

eo(n)  = ~Fi(?,c°)  (^o/h).  (5.56) 

As  n e must  conform  to  (5.38)  and  to  satisfy  this  requirement 
define 

(5.57) 

The  eddy  viscosity  defined  by  (5.56)  and  (5.57)  effectively  reverses 
the  behavior  inherent  in  (5.45);  in  this  case  the  limit  (5.42)  is 
approached  exponentially  quickly  but  the  decay  to  (5.38)  is  algebraic. 

In  data  comparisons  using  (5.56)  it  was  found  that  the  relatively  slow 
decay  of  (5.56),  for  large  n,  tended  to  degrade  the  quality  of  the  fits. 
Consequently  an  alternative  form  was  defined  according  to 


A 
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^ -L  ajQj(a„/n) 

J ^ 


Here  the  Q.  are  a set  of  functions  (to  be  selected)  which  decay  exponen- 

si 

tially  to  zero  as  n ■»■  0.  The  sum  of  terms  in  (5.58)  is  added  solely  to 

increase  the  rate  of  decay  to  the  limit  (5.38)  for  large  n and  the 

constants  a-  are  to  be  selected  for  this  purpose.  There  is  a wide 
%) 

variety  of  functions  Q.  which  will  suit  the  above  purpose  and  it  was 

sJ 

arbitrarily  decided  to  take  solutions  of  the  equation 


y"  + 2xy'  - 2(j  - l)y  = 0, 


(5.59) 


for  the  Q..  These  solutions  are  parabolic  cylinder  functions  and  may 

si 

be  related  to  repeated  integrals  of  the  error  function.  They  represent 
a complete  set  of  functions  and  the  first  few  Q.  are  defined  by: 

Qi(x)  = eAi(  X - 1,  ] 


Q,(x)  = x(eAj(  X - 1)  + 


2x  -x' 


Q3(x)  = (2x^  + l)(eA^  X - 1)  + ^ e 


(5.60) 


Q4(x)  = (2x^  + 3x)(eA|(  X - 1)  + — (x^  + l)e”^  , 

/tt 


Q,(x)  = (4x‘'  + 12x2  + 3)(eA^  x - 1)  + (2x'  + 5x)e■^^ 


In  Table  3,  the  coefficients  associated  with  retaining  an  increasing 
number  of  terms  in  the  sum  in  (5.58)  is  given  as  well  as  the  rate  of 
decay  to  the  limit  (5.38).  In  Figure  2 the  functions  C|j(n)  and  ^^(n) 
are  plotted.  It  is  clear  that  by  taking  additional  terms  in  (5.58) 
the  rate  of  decay  to  the  limit  (5.38)  may  be  accelerated;  however  the 
quality  of  the  fits  obtained  using  ec(n)  did  not  improve  substantially 
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Table  3 

Coefficients  Associated  with  the  Eddy  Viscosity 
Formula  (5.58)  for  Various  Values  of  N 


N 

0,  1,  2 

3 

4 

5 

a^K/K 

0.886227 

1.032408 

1.163729 

1.282963 

0 

1 

3.476918 

9.702465 

0 

3.544908 

18.488031 

74.728587 

^3 

0 

1 

10.430755 

68.268486 

a. 

0 

0 

3.081339 

43.304199 

^5 

0 

0 

0 

4.351232 

Rate  of  decay 
to  (5.38) 

oin^y 

0 

1 

0 

1 

0 

1 

cn 

over  that  obtained  using  e^(n).  Consequently,  formulae  for  Cj^(n)  for 
N > 5 were  not  considered. 

The  data  comparisons  using  c^{r])  are  given  in  Appendix  D for  the 
same  data  sets  considered  in  Appendix  C;  the  fits  described  in  Appendix 
D are  one-parameter  optimizations  on  the  parameter  u^  and  may  be  com- 
pared with  the  corresponding  cases  in  Appendix  C.  Overall  for  zero  and 
favorable  pressure  gradients,  use  of  £5(11)  produces  slightly  better 
results  than  use  of  equation  (5.45).  For  the  mild  adverse  pressure 
gradient  equilibrium  flow  of  Andersen,  et  al.  (1972),  there  is  a 
significant  improvement  in  the  quality  of  the  fits.  Thus  it  appears 
that  maintaining  the  strictly  linear  behavior  in  e(n)  near  n = 0 is 


( 


3 


I 


i 

i 

I 

1 


worthwhile. 
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5.6  Summary 

It  has  been  demonstrated  in  this  section  that  the  one-parameter 
composite  profile  consisting  of  the  unsteady  wall  layer  model  and  an 
outer  region  similarity  profile  gives  a good  representation  of  equi- 
librium flow  data.  In  zero  and  favorable  pressure  gradient  flows,  the 
quality  of  the  fits  improves  downstream,  as  expected  from  theoretical 
considerations.  In  the  adverse  pressure  gradient  flow  this  trend  is 
not  observed.  By  varying  the  eddy  viscosity  parameter  K as  well  as 
u^,  the  quality  of  the  fits  (at  data  stations  where  the  flow  is  not 
self-similar)  may  be  significantly  improved.  Lastly,  the  eddy  viscosity 
formula  used  in  §5.5  (wherein  an  essentially  linear  variation  is  imposed 
near  n = 0)  appears  to  be  somewhat  preferable  to  the  form  used  in  §5.4. 

One  final  point  is  worthy  of  note  here.  In  carrying  out  the  pro- 
file fits,  the  outer  region  length  scale  was  selected  as  6*/u*  according 
to  equation  (5.44).  The  main  reason  for  this  choice  is  that  the  eddy 
viscosity  formulae  (5.45)  and  (5.58)  assume  a simple  form  with 
f^i(x,“)  = -1.  However  an  equally  valid  choice  for  Ag  would  be  the 
boundary-layer  thickness  which  is  also  0(u*).  Define  the  boundary-layer 
thickness  as  that  ordinate  where  the  velocity  profile  achieves  some 
fixed  and  consistent  fraction  0 of  the  mainstream  value;  typical  defini- 
tions of  6 are  0.99  or  0.995.  From  equation  (5.3) 


u* 


9F, 

9n 


= -(1  - 0). 


(5.61) 


With  Aj  selected  according  to 
Ao  = 


6 


(5.62) 
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equations  (5,18)  and  (5.40)  become,  respectively. 


i = 1 iogjRej  u*U_(x)}  + C,  - C,, 


An  analysis  similar  to  that  leading  to  (5.50)  leads  to, 

+ 2B)n*^+23Uj  = 0 


(5.63) 

(5.64) 


(5.65) 


u*^62  dn*  dn*/  ' dr? 

where  q*  = y/6  and  8 is  the  same  as  in  equation  (5.50).  The  fitting 
procedure  described  in  §5.4  and  §5.5  was  also  carried  out  using  equa- 
tion (5.65).  To  initiate  the  procedure  a value  of  6 was  first  guessed 
and  a solution  of  (5.65)  computed;  the  value  of  6 was  refined  in  an 
iterative  procedure  subject  to  the  constraint  (5.61).  In  this  way 
results  identical  to  those  described  in  §5.4  and  §5.5  were  obtained. 
Clearly  the  two  approaches  are  equivalent  since  (5.50)  may  be  obtained 
from  (5.65)  by  a transformation  of  n*. 


kk^i 
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6.  SUMMARY  AND  CONCLUSIONS 

The  stated  objective  of  this  report  was  to  present  data  compari- 
sons for  four  models  for  the  velocity  profile  in  a two-dimensional, 
time-mean,  constant  property  turbulent  boundary  layer.  The  primary 
motivation  for  this  was  to  ascertain  how  well  each  model  represents 
experimental  data  and  to  offer  a critical  evaluation  of  each  model  on 
this  basis.  A secondary  goal  was  to  search  for  any  obvious  trends  in 
the  profile  parameters  which  might  prove  useful  in  a boundary  layer 
prediction  method. 

In  section  3,  the  Van  Driest  (1956)  and  extended  Van  Driest 
models  for  the  inner  layer  were  compared  with  the  more  recent  unsteady 
wall  layer  model  (Walker  & Abbott,  1977,  Walker  & Scharnhorst,  1977). 

To  provide  a direct  test  of  the  two  models  only  inner  region  data 
was  considered  and  the  results  are  described  in  Appendix  A.  The  data 
comparisons  and  the  analysis  in  section  3 show  the  failure  of  the  ex- 
tended Van  Driest  model  in  flows  with  pressure  gradient  and  for  this 
reason,  it  was  recommended  that  this  model  be  abandoned.  On  the  other 
hand,  the  unsteady  wall  layer  model  represents  the  data  very  well  even 
in  severe  pressure  gradients  and  generally  gives  a much  closer  repre- 
sentation of  the  data;  of  course  this  does  not  prove  the  validity  of 
the  assumptions  behind  the  unsteady  wall  layer  model  (Walker  & Abbott, 

1977  and  Walker  & Scharnhorst,  1977)  but  does  provide  a good  degree 
of  confidence  in  the  resulting  profile. 

From  a practical  standpoint  a number  of  points  in  connection 
with  the  unsteady  wall  layer  model  deserve  brief  consideration.  First  | 
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the  unsteady  wall  layer  model  contains  the  parameters  T and  which 
are  physically  observeable  quantities  of  the  flow;  this  is  in  con- 
trast to  the  parameter  A^,  in  the  Van  Driest  (1956)  and  extended  Van 
Driest  models,  which  is  simply  an  adjustable  constant  with  no  physical 
significance.  At  present  there  are  substantial  experimental  problems 
in  measuring  the  quantities  T and  u^  directly  but  it  may  be  expected 
that  these  difficulties  will  be  resolved  at  some  future  date.  In  the 
second  place,  the  theoretical  basis  of  the  unsteady  wall  layer  model 
is  somewhat  more  satisfactory  than  for  the  Van  Driest  (1956)  model. 

At  the  same  time,  there  are  a number  of  unresolved  features  surround- 
ing the  theoretical  development  of  the  unsteady  wall  layer  model 
which  are  discussed  in  Walker  & Scharnhorst  (1977).  Further  work  in 
this  area  is  clearly  needed  with  particular  reference  to  obtaining 
a good  theoretical  understanding  of  the  dynamics  of  the  outer  layer 
and  the  observed  inviscid-viscous  interactions  with  the  inner  layer. 

For  the  present,  the  profile  (3.16)  must  be  regarded  as  only  a 
useful  model.  Finally  equation  (3.16)  appears  more  complicated  than 
the  Van  Driest  (1956)  model;  however  the  computation  time  necessary  to 
perform  the  quadrature  in  (3.16)  accurately  is  comparable  to  that 
necessary  to  evaluate  (3.4).  In  addition  a program  to  evaluate  (3.16) 
is  available  from  the  authors  upon  request.  In  summary,  for  the  above 
reasons  the  unsteady  wall  layer  model  is  believed  to  be  a preferred 
and  useful  formula  for  the  mean  profile  in  the  wall  region. 

Additional  verification  of  the  unsteady  wall  layer  model  was  con- 
sidered in  sections  4 and  5,  where  (3.16)  was  coupled  with  Coles'  (1956) 
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"law  of  the  wake"  and  an  outer  region  similarity  profile  respectively, 
to  form  a composite  profile  for  the  entire  boundary  layer.  We  first 
consider  the  results  of  section  4 and  discuss  the  "law  of  the  wake". 

The  three  parameter  fits  using  (4.4),  have  been  shown  to  yield 
a close  representation  of  data  under  a wide  range  of  experimental 
conditions.  In  this  context  equation  (4.4)  may  be  considered  to  be 
very  useful  and  the  results  described  in  Appendix  B provide  addition- 
al confidence  in  the  unsteady  wall  layer  model.  Moreover,  the  data 
fits  also  suggest  that  assuming  that  has  a universal  value  of  5.0 
(or  some  other  value)  is  not  justified. 

However,  any  conclusions  based  on  use  of  the  "law  of  the  wake" 
must  be  regarded  with  some  caution  since  a number  of  legitimate  ob- 
jections may  be  raised  to  Coles  (1956)  model.  The  first  of  these  is 
that  the  model  contains  the  ill-defined  boundary  layer  thickness  6; 
in  data  comparisons  in  this  report  as  v/ell  as  in  Coles  & Hirst  (1969) 
this  quantity  must  be  treated  as  a parameter  to  obtain  good  repre- 
sentations of  the  data.  However,  the  optimized  values  of  6 do  not 
appear  to  correspond  to  a consistent  definition  of  the  boundary  layer 
thickness  from  data  set  to  data  set.  A second  theoretical  deficiency 
is  that  the  slope  and  all  higher  derivatives  of  (4.4)  are  not  contin- 
uous at  the  boundary  layer  dege.  The  third  and  most  serious  difficulty 
with  the  "law  of  the  wake"  is  the  notion  that  the  velocity  profile  in 
the  outer  layer  may  be  written  as  a linear  combination  of  a logarithm 
and  a wake  function;  this  deficiency  has  been  recently  discussed  by 
Fendell  (1972)  and  is  now  recognized  to  be  a flaw  in  the  original 
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formulation.  In  zero  and  adverse  pressure  gradient  flows  the  wake 
strength  parameter  assumes  relatively  large  values  which  tends  to 
mask  this  basic  problem.  However,  in  favorable  pressure  gradient 
flows,  the  wake  strength  parameter  becomes  small  and  it  is  in  these 
situations  where  the  improper  formulation  is  most  evident. 

Although  there  are  a number  of  well  known  disadvantages  inherent 
in  momentum  integral  prediction  methods,  the  composite  profile  (4.4) 
could  in  principle  be  used  for  this  purpose.  One  of  the  better  in- 
tegral methods  is  due  to  Alber  (1968)  which  has  been  recently  described 
by  White  (1974,  p.289).  This  method  has  achieved  a fair  degree  of 
success  and  is  based  on  Coles  (1956)  model.  Alber  (1968)  ignores  the 
inner  layer  and  assumes  a universal  value  for  C^. ; the  "law  of  the  wake" 
thus  provides  a two  parameter  expression  for  the  velocity  profile  in 
the  outer  layer.  However,  it  is  clear  from  Appendix  B that  a univer- 
sal value  of  C^.  will  not  provide  a close  representation  of  profile  data 
in  every  situation.  The  effect  of  the  inner  layer  could  be  included 
through  use  of  (4.4);  however,  (4.4)  is  a three  parameter  expression 
and  this  means  that  at  least  two  experimental  correlations  of  the 
parameters  are  required.  As  White  (1974,  p.482)  has  shown  the  wake 
strength  parameter  n can  be  correlated  as  a function  of  the  Clauser 
parameter  (defined  in  equation,  5.30)  and  such  a correlation  could 
be  obtained  using  the  results  of  Appendix  B.  However,  the  form  of  the 
required  second  correlation  relating  the  profile  parameters  is  not 
evident  at  this  time.  It  is  also  worthwhile  to  remark  that  there  is 
serious  question  about  the  validity  of  the  type  of  approach  used  by 
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Alber  (1968).  The  so-called  mechanical  energy  equation  or  its  inte- 
grated form  (containing  the  dissipation  integral)  is  not  a relation 
which  is  independent  from  the  momentum  equation.  For  these  reasons 
and  because  of  the  difficulties  discussed  in  connection  with  the 
"law  of  the  wake"  it  is  not  clear  that  this  momentum  integral  ap- 
proach can  be  pursued  with  profit. 

In  summary,  we  regard  the  "law  of  the  wake"  as  an  interesting 
correlation;  because  it  represents  outer  region  data  well,  it  is  a 
useful  expression  to  examine  the  validity  of  the  unsteady  wall  layer 
model.  Given  profile  data,  equation  (4.4)  may  be  used,  through  a 
three  parameter  optimization,  to  closely  represent  the  data  at  that 
station.  However,  because  of  various  deficiencies  in  the  "law  of 
the  wake"  any  attempt  to  use  the  formula  in  a prediction  method 
should  be  exercised  with  caution. 

In  sect"‘.'n  5,  a similarity  profile  was  used  for  the  outer  region. 
The  theoretical  foundation  for  this  profile  is  somewhat  firmer  than 
for  the  "law  of  the  wake"  and  good  representations  of  equilibrium  flow 
data  were  obtained  with  a one  parameter  (u^)  optimization  procedure. 

As  expected  the  profiles  in  zero  and  favorable  pressure  gradients  are 
successively  represented  better  at  subsequent  stations  downstream;  the 
adverse  pressure  gradient  flow  does  not  exhibit  this  trend  suggesting 
that  in  this  environment  similarity  is  not  being  approached.  It  was 
also  shown  that  in  situations  where  the  flow  is  not  similar;  a closer 
representation  of  the  data  could  be  achieved  by  varying  the  eddy  vis- 
cosity constant  K as  well  as  u^ . This  is  important  for  a differential 
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prediction  method,  where  the  data  at  some  initial  station  must  be  re- 
presented in  order  to  initiate  a numerical  marching  scheme  downstream. 
A prediction  method  based  on  solving  equation  (5.6)  numerically  in  the 
outer  layer  whilst  matching  to  the  inner  profile  (3.16)  is  currently 
under  consideration  and  will  be  reported  in  due  course. 
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Introduction 

Data  comparisons  between  the  Van  Driest  inner  region  model  given 
by  equation  (3.4)  and  the  unsteady  wall  layer  model  given  by  equation 
(3.16)  are  discussed  in  this  appendix.  The  data  sets  are  grouped 
according  to  pressure  gradient  and  a few  comments  ^i^e  presented  below. 

A.l  Zero  Pressure  Gradient  Data 

The  Van  Driest  and  unsteady  wall  layer  fits  to  the  constant  pressure 
data  of  Andersen,  et  al.  (1972)  are  detailed  in  Table  A.l  and  Figure  A.l 
and  Table  A. 2 and  Figure  A. 2,  respectively.  Visually  there  appears  to  be 
little  to  choose  between  the  two  models.  However  from  the  tables  e = 
0.0072  for  the  Van  Driest  model  while  e = 0.0073  for  the  unsteady  wall 
layer  profile;  here  e is  the  average  root-mean-square  error  over  the  nine 
data  stations.  Thus  on  average  the  Van  Driest  profile  appears  to  fit  the 
profiles  very  slightly  better;  however  the  difference  is  not  significant 
and  no  conclusions  can  be  drawn.  For  example,  both  profiles  do  not  fit 
the  sixth  data  station  very  well  and  if  this  data  station  is  omitted  e 
for  the  unsteady  wall  layer  model  is  slightly  lower  (e  = 0.0065)  than 
that  for  the  Van  Driest  model  (e  = 0.0066).  It  should  be  noted  that  in 
Table  A.l  the  optimum  value  for  A^  varies  substantially  about  the  value 
of  A^  = 26.0  suggested  by  Van  Driest  (1956).  It  may  also  be  observed 
that  the  optimum  value  of  u^  is  appreciably  different  from  the  quoted 
experimental  value  for  the  Van  Driest  model  than  for  the  unsteady  wall 
layer  model.  However  because  of  the  uncertainties  that  have  already 
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been  discussed  in  connection  with  the  methods  of  determining  from 
experimental  data,  no  definite  conclusions  can  be  made  here. 

A. 2 Mild  Adverse  Pressure  Gradient  Data 

The  fits  of  the  two  inner  region  models  to  a mildly  adverse  gradi- 
ent flow  are  shown  in  Figures  A. 3 and  A. 4 and  the  numerical  values  of 
the  parameters  are  given  in  Tables  A. 3 and  A. 4.  A close  inspection  of 
Figures  A. 3 and  A. 4 reveals  that  the  unsteady  wall  layer  model  appears 
to  give  a closer  representation  of  the  data  than  the  extended  Van  Driest 
profile  and  this  is  borne  out  from  the  tables  where  e = 0.0058  for  the 
unsteady  wall  layer  profile  and  e = 0.0076  for  the  extended  Van  Driest 
model.  The  extended  Van  Driest  profile  only  gives  a better  fit  on  the 
fifth  data  station  although  it  is  evident  from  the  tables  that  both 
profiles  encounter  difficulty  for  this  station.  For  this  data  set  the 
unsteady  wall  layer  model  is  superior  and  the  difference  is  significant. 

A. 3 Strong  Adverse  Pressure  Gradient  Data 

The  results  of  the  two  profile  fits  to  inner  region  data  for 
strongly  adverse  gradient  flows  are  shown  in  Figures  A. 5 to  A. 8 and 
Tables  A. 5 to  A. 8;  these  include  the  third  Andersen,  et  al . (19T2)  flow 
and  Fraser  Flow  A (Coles  & Hirst,  1969). 

The  data  fits  for  the  Andersen,  et  al.  (1972)  flow  are  given  in 
Figures  A. 5 and  A. 6;  this  flow  has  been  referred  to  as  a strong  adverse 
pressure  gradient  flow  by  the  originators  since  the  values  of  p^  are 
about  twice  as  large  as  for  their  second  flow  (discussed  in  A. 2).  On 
average  the  unsteady  wall  layer  model  fits  the  data  (e.  = 0.0062)  better 
than  the  extended  Van  Driest  model  (c  = 0.0070).  The  extended  Van  Driest 
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model  gives  a better  fit  on  only  the  fourth  and  sixth  data  stations. 

The  data  comparisons  for  Fraser  Flow  A are  given  in  Figures  A. 7 
and  A. 8;  for  this  flow  the  values  of  p’*'  are  larger  than  in  the  previous 
case  and  in  fact  this  flow  is  proceeding  toward  separation.  At  the 
first  two  stations  the  value  of  p'*'  is  small;  both  profile  models  give  a 
good  representation  of  the  data  although  e is  smaller  for  the  unsteady 
wall  layer  model.  At  subsequent  data  stations  downstream  as  p'*' 
increases,  the  extended  Van  Driest  model  begins  to  deviate  drastically 
from  the  data.  Recall  that  the  result  of  the  optimization  procedure  is 
the  optimum  performance  for  a model  (according  to  the  fitting  criterion 
(1.5)).  Clearly  there  is  no  correspondence  between  the  extended  Van 
Driest  model  and  these  data;  the  algebraic  behavior  of  the  profile 
(discussed  in  §3.1)  is  clearly  evident  in  Figure  A. 7 and  the  result  of 
this  is  that  the  data  can  only  skew  through  the  data  points.  On  the 
other  hand,  the  unsteady  wall  layer  model  conforms  to  the  latter  data 
stations  quite  well  as  is  apparent  in  Figure  A. 8;  there  is  no  drastic 
deviation  from  the  data  as  in  Figure  A. 7.  The  c for  the  extended  Van 
Driest  model  is  almost  three  times  as  large  as  that  for  the  unsteady 
wall  layer  model  (0.0181  versus  0.0066). 
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APPENDIX  B 

DATA  COMPARISONS  FOR  FULL  PROFILE  DATA  USING  A 
COMPOSITE  PROFILE  CONSISTING  OF  THE  UNSTEADY 
WALL  LAYER  MODEL  AND  THE  LAW  OF  THE  WAKE 
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Introduction 

In  this  appendix,  full  profile  data  are  considered  and  equation 
(4.4)  is  fit  by  optimization  of  the  three  parameters  u^,  6 and  C^.  The 
results  are  presented  in  graphical  and  tabular  form.  A few  comments 
about  each  data  set  follow. 

1300  - Ludwieg  and  Tillman 

The  fits  to  the  accelerating  flow  of  Ludwieg  and  Tillman  are 
described  in  Table  B.l  and  Figure  B.l.  The  overall  representation  of 
the  data  by  the  composite  model  is  very  good  with  the  worst  fits 
occurring  at  stations  2,  3 and  9.  The  negative  and  near  zero  values 
of  the  wake  strength  parameter  n influence  the  quality  of  the  fits 
adversely  at  stations  2 and  3.  As  stated  by  Coles  and  Hirst  (1969)  the 
data  at  station  9 are  distorted  and  may  not  be  reliable.  The  flow  is 
nearly  in  equilibrium  so  that  II  should  be  approximately  constant;  a 
fair  estimate  is  about  0.2  over  the  middle  stations.  Note  the  trend 
for  and  S to  increase  downstream. 

2100  - Schubauer  and  Klebanoff 

The  comparisons  to  the  data  of  Schubauer  and  Klebanoff,  taken  on  a 
large  airfoil-like  body,  are  given  in  Table  B.2  and  Figure  B.2.  Only 
the  last  fifteen  stations  are  considered  here  since  they  represent  that 
part  of  the  flow  where  the  pressure  gradient  is  most  severe  as  the  flow 
proceeds  toward  separation.  Examination  of  the  root-mean-square  error 
in  Table  B.2  suggests  that  the  fits  are  only  fair.  This  is  believed 
due  to  two  reasons.  In  some  cases  there  are  not  enough  data  points  to 


clearly  define  the  boundary-layer  edge.  Second,  some  of  the  data  points 
near  the  wall  show  an  erratic  behavior  at  some  stations,  which  adversely 
influences  the  fitting  procedure  and  the  root-mean-square  error.  Indeed 
the  quality  of  some  of  the  data  is  suspect  as  stated  in  Coles  and  Hirst 
(1969).  However,  note  the  trend  of  the  profile  parameters;  as  the 
pressure  gradient  parameter  increases  both  6 and  C^.  decrease. 

2300  - Clauser  2 

The  fits  to  the  moderately  adverse  pressure  gradient  flow  of 
Clauser  are  described  in  Table  B.3  and  Figure  B.3.  The  fits  to  this 
equilibrium  flow  are  only  fair.  Part  of  this  may  be  due  to  a lack  of 
inner  region  points  near  the  wall.  Coles  and  Hirst  (1969)  also  remark 
that  this  flow  is  slightly  three  dimensional.  Note  that  the  IT  parameter 
is  roughly  constant  at  about  4. 

2500  - Bradshaw  Equilibrium,  a = -0.15 

The  fits  to  this  mildly  adverse  gradient,  equilibrium  flow  are 
described  in  Table  B.4  and  Figure  B.4.  The  fits  are  fair.  Note  that 
n is  roughly  equal  to  about  1.3  for  all  four  profiles,  while  both  6 and 
C^.  decrease.  It  appears,  upon  examination  of  Figure  B.4,  that  the  slope 
of  a line  through  the  data  points  in  the  logarithmic  region  is  somewhat 
different  from  0.41. 

2600  - Bradshaw  Equilibrium,  a = -0.255 

The  fits  to  this  moderately  adverse  gradient,  equilibrium  flow  are 
described  in  Table  B.5  and  Figure  B.5.  The  fits  are  not  very  good. 

Note  the  tendency  of  the  profile  to  skew  through  the  data  in  the  "law 
of  the  wall"  region  suggesting  a k different  from  0.41.  There  are  some 
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unresolved  difficulties  discussed  by  Coles  and  Hirst  (1969)  with  this 
data  set  and  also  with  flow  2500. 

2700  - Herring  and  Norbury,  g = -0.35 

The  fits  to  this  mildly  favorable  gradient  equilibrium  flow  are 
described  in  Table  B.6  and  Figure  B.6.  Overall  the  fits  are  good. 

Note,  however,  that  as  with  flow  1300,  a near  zero  or  negative  n value 
does  not  produce  a good  fit  to  the  data.  This  situation  occurs  at 
stations  3 and  4 where  the  quality  of  the  fit  is  significantly  degraded. 

The  inner  region  data  are  fit  quite  well  though. 

2800  - Herring  and  Norbury,  B = -0.53 

The  fits  to  this  strongly  favorable  gradient  equilibrium  flow  are 
described  in  Table  B.7  and  Figure  B.7.  The  fits  are  all  excellent  even 
in  the  case  of  small  and  negative  II  values.  The  wake  component  of  the 
outer  profile  diminished  as  n decreases  and  the  fact  that  n becomes  small 
and  negative  at  the  latter  data  stations  may  be  indicative  of  relaminari- 
zation.  Indeed,  Coles  and  Hirst  (1969)  suggest  this  may  have  occurred  | 

in  this  flow.  In  spite  of  this,  the  composite  profile  represents  the 
data  very  well. 

2900  - Perry 

The  data  comparisons  for  this  diverging  channel  flow  are  described 
in  Table  B.8  and  Figure  B.8.  Overall  the  fits  to  this  adverse  gradient 
diffuser-like  flow  are  very  good.  The  worst  fits  occur  at  the  third  and  ■ 

fourth  stations  where  it  is  apparent  that  the  data  are  incomplete  near 
the  boundary-layer  edge.  The  composite  profile  represents  the  data  very 
well  as  p"*"  increases  even  through  the  last  two  stations  where  flow 

^ 


separation  is  imminent;  observe  the  trend  for  S and  C^.  to  decrease  down- 
stream. Also  it  may  be  observed  from  Figure  B.8  that  there  is  a slight 
tendency  for  the  profile  to  skew  through  the  logarithmic  zone  which 
suggests  a value  of  k different  from  0.41  may  be  appropriate. 

3300  - Bradshaw  Layer  C 

The  fits  for  this  case  are  given  in  Table  B.9  and  Figure  B.9.  The 
flow  is  initially  at  constant  pressure  and  is  developing  into  equilibrium 
in  a moderate  adverse  pressure  gradient.  The  fits  for  this  case  are 
excellent  with  the  smallest  e recorded  in  this  report  achieved  at  station 
2.  The  n parameter  increases  at  the  downstream  stations  and  appears  to 
be  approaching  an  approximately  constant  value  as  the  flow  approaches 
equilibrium.  A similar  trend  for  S and  may  be  observed  in  Table  B.9. 

3500  - Newman  Airfoil  - Series  2 

The  fits  to  this  airfoil  boundary-layer  flow  proceeding  to  separa- 
tion are  described  in  Table  B.IO  and  Figure  B.IO.  In  general  the  fits 
are  quite  good.  As  p'*'  increases  downstream  S and  C^-  decrease  and  n 
increases.  Note  that  at  the  last  station,  which  is  believed  to  be  very 
close  to  separation,  the  logarithmic  region  has  apparently  disappeared. 

3600  - Moses  1 

The  fits  to  this  flow  are  described  in  Table  B.ll  and  Figure  B.ll. 
This  case  is  the  boundary  layer  on  a cylinder  in  an  axially  symmetric 
flow.  The  rate  of  pressure  increase  is  initially  slow  but  then  quite 
abrupt.  In  general  the  quality  of  the  fits  is  good.  The  sparseness 
of  the  inner  region  makes  it  difficult  to  judge  the  inner  model  fit; 
however  the  Coles'  model  fits  the  data  well.  Observe  that  is 


significantly  different  from  5.0  and  decreases  as  separation  is 
approached;  the  last  station  is  evidently  close  to  separation. 

3700  - Moses  2 

The  fits  for  this  case  are  given  in  Table  B.12  and  Figure  B.12. 

This  case  corresponds  to  the  boundary  layer  on  a cylinder  in  an  axially 
symmetric  flow;  the  rate  of  pressure  rise  is  moderate.  The  overall 
quality  of  the  fits  is  fair.  Again  there  is  a sparseness  of  inner  region 
points  which  makes  it  difficult  to  assess  the  wall  layer  model.  Note 
that  C^.  is  significantly  different  from  5.0  and  also  that  for  the  first 
few  stations,  the  analytical  profile  skews  slightly  through  the  logarith- 
mic region  data.  This  suggests  a value  of  k different  from  0.41. 

3800  - Moses  3 

The  fits  to  the  third  Moses  cylinder  flow  are  given  in  Table  B.13 
and  Figure  B.13.  For  this  case  there  is  a strong  pressure  rise.  The 
last  two  data  stations  are  apparently  in  a separated  flow  region  and 
were  omitted  from  consideration  for  this  reason.  The  fits  are  fair,  at 
best.  For  the  last  three  profiles  fitted,  the  analytical  curve  skews 
through  the  logarithmic  zone  suggesting  a value  of  k different  from  0.41 
may  be  appropriate. 

4000  - Moses  5 

The  fits  for  the  fifth  Moses  cylinder  flow  are  given  in  Table  B.14 
and  Figure  B.14.  In  this  case  there  is  an  initially  strong  pressure  rise 
followed  by  a relaxation  at  constant  pressure.  It  is  apparent  that  the 
fits  are  not  good  for  this  flow.  The  sparseness  of  the  inner  region 
data,  coupled  with  the  fact  that  the  data  points  that  are  in  the  inner 


region  appear  suspect,  makes  it  difficult  to  draw  any  conclusions  about 
the  fits.  In  actuality,  smaller  values  of  e were  obtained  for  this  data 
set  but  the  optimal  boundary-layer  thickness  6 took  on  unusually  large 
values.  The  fits  presented  here  are  those  which  appear  to  be  the  most 
reasonable. 

4400  - Schubauer  and  Spangenberg  A 

The  fits  to  this  strong  positive  pressure  gradient  flow  are  given 
in  Table  B.15  and  Figure  B.15.  It  is  evident  that  the  fits  are  not  good 
for  this  case.  The  same  problem  that  occurred  for  flow  4000  in  connec- 
tion with  6 was  encountered  in  this  case  as  well.  There  are  not  enough 
points  to  properly  define  the  inner  region  and  possible  errors  in  the 
data  points  near  the  wall  may  adversely  affect  the  fits.  Another  possible 
difficulty  here  is  the  truncation  of  the  profiles  disc.'ssed  by  Coles  and 
Hirst  (1969). 

4500  - Schubauer  and  Spangenberg  B 

The  fits  to  this  moderate  positive  pressure  gradient  flow  are  given 
in  Table  B.16  and  Figure  B.16.  The  fits  for  this  case  are  fair.  Some 
of  the  data  near  the  wall  appears  somewhat  suspect  which  tends  to  degrade 
the  quality  of  the  fit.  The  composite  model  represents  the  data  reason- 
ably well  up  to  the  last  station  which  is  close  to  separation.  The  trend 
for  S and  C^.  to  decrease  downstream  may  be  observed  in  Table  B.16. 

5000  - Fraser  A 

The  comparisons  for  this  conical  diffuser  flow  are  given  in  Table 
B.17  and  Figure  B.17.  Overall  the  fits  to  this  adverse  gradient  flow 
are  quite  good.  The  largest  values  of  e occur  at  stations  4,  5 and  10. 
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The  mainstream  velocity  gradient  is  changing  rapidly  at  the  fourth 
and  fifth  stations;  consequently  the  quoted  values  of  the  mainstream 
velocity  gradient  (in  Coles  & Hirst,  1969)  may  be  somewhat  in  error. 

The  major  portion  of  the  error  at  the  tenth  station  derives  from  the 
first  two  inner  region  points  and  the  last  few  outer  layer  points.  It 
is  interesting  to  observe  the  gradual  disappearance  of  the  logarithmic 
region  as  the  flow  approaches  separation  shortly  after  the  last  station. 

5100  - Fraser  B 

The  fits  to  Fraser's  B flow  are  similar  to  those  of  flow  A dis- 
cussed above  as  evident  from  Table  B.18  and  Figure  B.18.  Again  the 
worst  fits  occur  at  the  fourth  and  fifth  stations  where  the  velocity 
gradient  is  changing  rapidly.  The  pressure  gradient  parameter,  p"*", 
continuously  increases  and  eventually  becomes  0(1)  at  the  last  station. 
This  profile  exhibits  no  logarithmic  region  and  is  essentially  all 
wake.  The  composite  profile  fits  it  fairly  well  though. 

5200  - Stratford  5 

The  fits  to  this  flow  (in  a channel  with  an  abrupt  onset  of 
moderate  positive  pressure  gradient)  are  described  in  Table  B.19  and 
Figure  B.19.  At  first  the  results  appear  disturbing,  especially  for 
the  second  and  third  stations.  However,  Coles  and  Hirst  (1969)  remark 
that  the  flow  at  these  stations  cannot  be  made  to  conform  to  the  law 
of  the  wall,  the  profiles  being  essentially  inviscid  in  character.  If 
these  two  profiles  are  discarded  the  fits  are  fair.  Note  that  these 
data  also  suggest  a value  of  k different  from  0.41  may  be  appropriate. 
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5300  - Stratford  6 

This  flow  is  a boundary-layer  flow  in  a channel  with  the  abrupt 
onset  of  a severe  positive  pressure  gradient.  The  fits  are  described 
in  Table  B.20  and  Figure  B.20  and  are  similar  to  those  obtained  for 
flow  5200.  Again  the  second  and  third  profiles  are  essentially  inviscid 
in  character.  The  raggedness  of  the  data  at  the  sixth  and  seventh 
stations  contribute  greatly  to  the  root-mean-square  error. 

7000  - Samuel  and  Joubert  (1974) 

The  fits  for  this  flow  in  an  increasingly  adverse  pressure  gradient 
are  described  in  Table  B.21  and  Figure  B.21.  The  fits  are  quite  good 
with  the  largest  root-mean-square  error  obtained  at  the  last  two  stations 
where  the  flow  is  undergoing  relaxation.  The  analytical  profile  skews 
slightly  through  the  logarithmic  region  suggesting  that  a value  of  ic 
different  from  0.41  may  be  appropriate. 

8100  - Andersen,  et  al.  (1972) 

The  excellent  fits  to  this  constant  pressure  flow  are  described 
in  Table  B.22  and  Figure  B.22.  The  large  number  of  inner  region  points 
provide  a good  definition  of  the  inner  region.  It  may  be  observed  that 
there  is  negligible  skewing  of  the  profile  through  the  logarithmic  zone 
and  that  all  but  the  first  two  profiles  appear  close  to  being  similar. 
Both  S and  IT  are  nearly  constant  for  the  last  seven  stations. 

8200  - Andersen,  et  al . (1972) 

The  fits  to  the  mildly  adverse  gradient  flow  of  Andersen,  et  al. 
(1972)  are  described  in  Table  B.23  and  Figure  B.23.  The  fits  are 
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excellent  although  again  the  analytical  profile  skews  slightly 
through  the  logarithmic  region.  The  fact  that  this  is  a nearly  equi- 
librium flow  is  indicated  by  the  approximately  constant  values  of 
S,  C^.  and  II. 

8300  - Andersen,  et  al.  (1972) 

The  fits  for  the  strongly  adverse  gradient  flow  of  Andersen, 
et  al.  (1972)  are  given  in  Table  B.24  and  Figure  B.24.  Again  the 
fits  are  excellent  but  a slight  skewing  of  the  profile  in  the 
logarithmic  zone  may  be  observed.  The  flow  is  near  equilibrium  as 
suggested  by  the  essentially  constant  values  of  II  and  S. 
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Figure  B.5.  Bradshaw  moderately  adverse  gradient  equilibrium 

flow,  a = -0.255.  (Note  shifted  origins  on  the  LT*" 
axis) . 
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Figure  B.16*  (continued). 
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APPENDIX  C 

DATA  COMPARISONS  USING  THE  COMPOSITE  PROFILE 
CONSISTING  OF  THE  UNSTEADY  WALL  LAYER  MODEL 
AND  THE  SIMILARITY  PROFILE. 

EDDY  VISCOSITY:  e = e^(n),  equation  (5.45). 
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Introduction 

The  data  comparisons  of  the  composite  profile  consisting  of  the 
unsteady  wall  layer  model  and  the  similarity  profile  (the  solution  of 
equation  (5.50))  are  described  in  this  appendix.  The  eddy  viscosity 
formula  used  here  is  given  by  equation  (5.45).  In  all  cases  the  values 
of  6*  used  in  the  fitting  were  taken  directly  from  experiment;  since 
the  displacement  thickness  is  an  integrated  quantity  obtained  from 
profile  measurements  the  experimental  percentage  error  inherent  in  this 
quantity  is  expected  to  be  small. 

The  composite  profile  (5.51)  contains  either  implicitly  or  explicitly 
the  three  independent  parameters  u^,  tc,  and  K.  With  k set  equal  to  0.41, 
as  in  the  'oles  (1956)  model  (Appendix  B)  and  K = 0.0168,  equation  (5.51) 
is  a one-parameter  family  of  profiles  with  u^  being  the  single  parameter. 
First  the  data  comparisons  will  be  carried  out  by  optimization  on  u^ 
alone;  in  the  latter  part  of  the  appendix  the  effect  of  varying  both  •: 
and  K will  be  considered. 

It  is  worthwhile  to  re-emphasize  that  the  outer  region  similarity 
profile  can  at  best  apply  to  equilibrium  flows;  moreover  in  view  of  the 
discussion  in  §5.1,  the  similarity  profile  should  give  a progressively 
better  representation  of  the  data  at  subsequent  stations  downstream  in 
constant  pressure  and  favorable  pressure  gradient  flows.  In  an  adverse 
pressure  gradient,  applied  over  a very  long  distance,  separation  is 
anticipated  at  some  downstream  location.  Consequently  a similarity 
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solution  is  not  expected  to  give  a good  representation  of  the  data  in 
an  adverse  pressure  gradient  flow  at  any  downstream  location. 

C.l  One-Parameter  Fits  (k  = 0.41,  K = 0.0168,  u^  varied) 

In  this  section  the  fits  obtained  by  varying  the  single  parameter 
u^  are  described.  Note  that  in  this  case  the  composite  profile  (5.51) 
contains  only  one  adjustable  parameter  whereas  the  profile  (4.4),  which 
uses  the  Coles  (1956)  model  and  which  was  treated  in  Appendix  B,  contains 
three  independent  parameters. 

1300  - Ludwieq  and  Tillman 

The  fits  of  the  composite  profile  (5.51)  to  this  favorable  pressure 
gradient  equilibrium  flow  are  described  in  Table  C.l  and  Figure  C.l.  In 
general  the  fits  are  good  but  overall  slightly  worse  than  those  obtained 
with  the  Coles  model  (Table  B.l),  particularly  at  the  upstream  stations. 
However,  the  trend  for  the  root-mean-square  error  to  decrease  at  subse- 
quent stations  downstream  is  encouraging;  in  fact,  at  the  downstream 
station  the  values  of  e become  smaller  than  those  obtained  with  the 
Coles  model  (Table  B.l). 

2700  - Herring  and  Norbury 

The  fits  to  this  mildly  favorable  equilibrium  flow  are  described  in 
Table  C.2  and  Figure  C.2.  The  quality  of  the  fits  is  good  and  the  data 
are  represented  better  on  average  than  with  the  Coles  (1956)  model 
(Table  B.6).  This  is  encouraging  in  view  of  the  fact  that  the  fits 
described  in  Table  B.6  are  three-parameter  fits. 
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8100  - Andersen,  et  al.  (1972) 

The  data  comparisons  for  this  constant  pressure  flow  are  described 
in  Table  C.3  and  Figure  C.3.  The  fits  are  not  as  good  as  obtained  with 
the  Coles  (1956)  model  in  Table  B.22.  However,  there  is  a trend  for  e 
to  decrease  downstream  which  is  encouraging;  moreover  there  are  two  less 
adjustable  parameters  here  than  in  Table  B.22. 

8200  - Andersen,  et  al.  (1972) 

The  fits  for  this  mildly  adverse  pressure  gradient  flow  are 
described  in  Table  C.4  and  Figure  C.4.  It  is  clear  that  these  one- 
parameter  fits  are  rather  poor  and  are  not  nearly  as  good  as  obtained 
in  Table  B.23  with  the  Coles  (1956)  model.  The  trend  of  e for  this  case 
is  revealing.  This  is  an  equilibrium  flow  but  in  an  adverse  pressure 
gradient  and  in  view  of  §5.1,  similarity  is  not  expected.  This  is  borne 
out  in  Table  C.4  where  the  fits  are  uniformly  poor;  that  is  to  say, 
there  is  no  trend  for  e to  decrease  significantly  at  the  downstream 
stations  which  was  observed  in  Tables  C.l  to  C.3. 

C.2  Two-Parameter  Fits  (k  = 0.41,  K varied,  u^  varied) 

In  this  section,  the  effect  of  varying  the  outer  region  eddy  vis- 
cosity constant  K as  well  as  u^  will  be  considered  for  the  zero  pressure 
gradient  flow  of  Andersen,  et  al.  (1972).  The  result  of  this  procedure 
is  described  in  Table  C.5  and  Figure  C.5.  The  motivation  for  this  pro- 
cedure is  two-fold.  In  the  first  place,  true  similarity  is  not  antici- 
pated in  this  flow  until  a large  distance  downstream;  thus  neglecting 
3Ui/9x  in  the  momentum  equation  is  not  strictly  justified.  Varying  K is 
equivalent  to  attempting  to  account  for  the  non-similar  behavior.  A 
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second  and  more  important  reason  concerns  the  ultimate  application  of 
this  theory  in  a prediction  method.  In  order  to  initiate  a prediction 
method  it  is  necessary  to  be  able  to  obtain  a good  analytical  representa- 
tion of  the  data  at  any  initial  station.  Allowing  K to  vary  provides  a 
convenient  mechanism  for  this  purpose. 

The  trend  for  K to  decrease  to  an  approximately  constant  value  at 
the  downstream  stations  may  be  observed  in  Table  C.5.  Moreover  it  is 
clear  from  Table  C.5  that  the  additional  parameter  significantly  improves 
the  quality  of  the  fit.  The  average  root-mean-square  error  is  reduced 
by  38%  from  the  one-parameter  fit  described  in  Table  C.3. 

C.3  Three-Parameter  Fits  (<  varied,  K varied,  u^  varied) 

In  this  section  all  three  parameters  were  varied  for  the  zero  pres- 
sure gradient  flow  of  Andersen,  et  al.  (1972).  The  motivation  for 
allowing  < to  vary  was  that  in  some  of  the  data  sets  a skewing  of  the 
profile  was  observed  through  the  logarithmic  zone.  This  tendency  may 
be  observed  by  close  inspection  of  Figure  C.5  and  suggests  a value  of  k 
different  from  0.41. 

The  result  of  this  optimization  procedure  is  described  in  Table 
C.6  and  Figure  C.6,  The  result  of  varying  the  third  parameter  is  to 
reduce  the  root-mean-square  error  from  Table  C.5  but  the  reduction  is 
not  as  substantial  as  was  obtained  between  Table  C.3  and  Table  C.5. 

The  fits  are  good  and  are  of  comparable  quality  to  that  obtained  with 
the  Coles  (1956)  model  in  Table  B.22.  The  primary  source  of  error  in 
this  case  originates  at  the  boundary-layer  edge  where  the  similarity 
profile  appears  to  decay  too  quickly  to  the  mainstream  value. 
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(Figure  C.2.  Herring  and  Norbury  mildly  favorable  pressure 

gradient  equilibrium  flow,  3=-0.35.  (Note 
shifted  origins  on  the  U"*"  eucis) . 
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TABLE  C.4:  ANDERSEN*  KAYS  AND  MOFFAT  - ILD  P,  G,  - 071571 

(ENGLISH  UNITS) 


Y+ 


Figure  C.4.  Mild  adverse  pressure  gradient  flow  of 

Andersen,  et  al.  (1972).  (Note  shifted  origins 
on  the  U'*'  axis)  . 
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APPENDIX  0 

DATA  COMPARISONS  USING  THE  COMPOSITE  PROFILE 
CONSISTING  OF  THE  UNSTEADY  WALL  LAYER  MODEL 
AND  THE  SIMILARITY  PROFILE. 

e = e5(n),  equation  (5.58). 
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Figure  D.l. 


Ludwieg  and  Tillman  accelerating  flow 
(Note  shifted  origins  on  the  U axis) 
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Figure  D.2.  Herring  and  Norbury  mildly  favorable 
pressure  gradient  equilibrium  flow, 
3=-0.35.  (Note  shifted  origins  on  the 

(;■*■  axis) . 
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